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1  OVERVIEW  OF  PROGRAM 

This  program  addresses  the  Navy  need  for  extended  firm  track  range  for  low  altitude  cruise 
missiles  through  the  integration  of  multiple  sensors.  Track-Before-Declare  (TBD)  techniques 
that  utilize  signal  features  are  proposed  for  the  synergistic  integration  of  an  Electronically 
Scanned  Array  (ESA)  radar  with  other  sensors  for  the  detection  of  weak  targets.  A  comparison 
of  the  performances  of  a  notional  multisensor  system  with  that  of  a  federated  system  of  sensors  is 
planned.  More  specifically,  the  integration  of  an  ESA  radar  and  Infrared  Search  and  Track 
(IRST)  sensor  for  a  shipboard  combat  system  was  proposed.  The  computer  simulation  models  of 
the  radars  and  IRST  sensor  will  include  the  effects  of  many  issues  such  as  finite  sensor 
resolution,  limitations  on  the  sensor  resources,  atmospheric  refraction,  sensor  pointing  errors, 
sea-surface  induced  multipath,  nonhomogeneous  clutter,  sea  clutter,  etc.  that  are  omitted  in  most 
of  the  legacy  simulations.  The  computer  simulation  models  will  be  utilized  to  develop  tracking 
benchmark  problems  for  broad  distribution.  These  benchmark  problems  will  serve  to  educate  the 
research  community  on  many  of  the  “real-world”  problems  that  are  faced  in  actual  tracking 
systems  and  the  integration  of  multiple  sensors. 

The  technical  objectives  of  the  project  include  the  following: 
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•  Development  and  demonstration  through  computer  simulation  of  algorithms  that  provide 
enhanced  detection  of  weak  targets  through  the  integration  of  multiple  sensors.  More 
specifically,  Track-Before-Declare  (TBD)  techniques  that  utilize  signal  features  are 
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proposed  for  the  synergistic  integration  of  an  Electronically  Scanned  Array  (ESA)  radar 
with  others  sensors  for  the  detection  of  weak  targets.  The  integration  of  the  ESA  radar 
with  the  Infrared  (IR)  sensor  are  to  be  considered.  The  detection  performance  of  a 
notional  multisensor  system  will  be  compared  with  that  of  the  federated  system  of 
sensors. 

•  Development  and  demonstration  through  computer  simulation  of  efficient  radar  resource 
allocation  techniques  that  maximizes  the  information  procured  by  the  multisensor  suite, 
while  accommodating  remote  cues  and/or  warnings  and  adapting  to  changes  in  the 
characteristics  of  the  targets  of  interest,  weather  conditions,  etc.  The  waveforms  and 
revisit  times  of  the  ESA  will  be  selected  to  maximize  the  information  procured  by  the 
multisensor  system.  The  detection  performance  of  a  notional  multisensor  system  with  the 
new  resource  allocation  techniques  will  be  compared  to  that  of  the  multisensor  system 
with  conventional  resource  allocation  and  that  of  the  federated  system  of  sensors. 

This  report  summaries  the  progress  and  accomplishments  made  during  October  1,  1999  and 
September  30,  2000,  which  is  the  second  year  of  a  three  year  program.  Since  only  a  small 
amount  of  funding  was  provided  for  this  program,  the  progress  toward  the  goals  of  this  program 
include  a  few  publications  and  other  related  activities.  Section  2  lists  the  papers  that  have  been 
published  or  submitted  for  publication,  and  Section  3  summarizes  related  activities  that  were 
supported  as  part  of  this  project. 

2  PUBLICATIONS 

The  following  four  papers  have  been  accomplished  through  funding  from  this  grant.  Copies  of 
the  papers  are  included  in  Enclosures  1  though  5. 

1 .  W.  D.  Blair  and  M.  Brandt-Pearce,  “Monopulse  DOA  Estimation  for  Two  Unresolved 
Rayleigh  Targets,”  IEEE  Trans,  on  Aerospace  and  Electronic  Systems,  to  appear  in  2001 . 

2.  G.C.  Brown,  W.D.  Blair,  and  D.A.  Diaz  “Track  Management  Technique  for  Electronically 
Scanned  Radars,”  in  Signal  and  Data  Processing  for  Small  Targets  2000,  Oliver  Drummond, 
Editor,  Proceedings  of  SPIE  Vol.  4048,  pp.  203-210  (2000) 

3.  W.  D.  Blair  and  M.  Brandt-Pearce,  “Effects  of  Diffuse  Multipath  on  the  Statistics  of 
Monopulse  Measurements,”  accepted  for  the  2001  IEEE  Radar  Conference,  Atlanta  GA, 

May  2001. 
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4.  W.  Wong  and  W.  D.  Blair,  “Steady-State  Tracking  With  LFM  Waveforms,”  IEEE 
Transactions  on  Aerospace  and  Electronic  Systems,  Vol.  36,  No.  2,  April  2000,  pp.  701-709. 

5.  W.  Wong  and  W.  D.  Blair,  “Steady-State  Tracking  With  LFM  Waveforms,”  Proceedings  of 
32nd  Southeastern  Symposium  on  Systems  Theory,  Tallahassee,  FL,  March  2000,  pp.  69-73. 

3  RELATED  ACTIVITIES 

The  follow  related  activities  have  been  supported  at  the  request  of  the  Office  of  Naval  Research. 

1.  Hosted  the  3rd  ONR/GTRI  Workshop  on  Target  Tracking  and  Sensor  Fusion,  May  17-18, 
2000  at  the  Cobb  Research  Facility  of  the  Georgia  Tech  Research  Institute.  The  workshop 
included  approximately  50  attendees  and  20  technical  presentations.  A  proceedings  of  the 
presentations  were  produced  as  a  bound  copy  and  on  CD  ROM. 

2.  Participated  as  Technical  Expert  in  Plenary  Panel  Discussion  entitled  “Fusion:  Vision  and 
Challenges,”  Fusion  2000:  3rd  International  Conference  on  Information  Fusion,  Paris  France, 
July  10-13  (see  Enclosure  6  for  slides  of  presentation). 

3.  Participated  in  several  meetings  with  personnel  from  the  United  Kingdom  on  a  Cooperative 
R&D  Program  on  Data  Fusion. 

>  February  21-22,  2000  in  NAVAIR  in  Lexington  Park,  MD:  Participated  in  discussions  on 
clutter  modeling,  data  fusion,  and  naval  weapons  systems  in  UK  and  USA  that  utilize 
data  fusion.  A  presentation  on  the  JCTN  Benchmark  was  given. 

>  September  14-15,  2000,  DERA  in  Portsmouth,  UK:  Participated  in  discussions  with  UK's 
DERA  personnel  on  the  UK's  R&D  programs  that  involve  data  fusion.  A  brief  entitled 
"Multiplatform/Multisensor  Tracking"  was  given  (see  Enclosure  7). 


Respectfully, 

William  Dale  Blair,  Senior  Research  Engineer 
Project  Director 
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ABSTRACT 

This  paper  provides  for  new  approaches  to  the  processing  of  unresolved  measurements  as 
two  DOA  measurements  for  tracking  closely-spaced  targets  rather  than  the  conventional  single 
DOA  measurement  of  the  centroid.  The  measurements  of  the  two  closely-spaced  targets  are 
merged  when  the  target  echoes  are  not  resolved  in  angle,  range,  or  radial  velocity  (i.e.,  Doppler 
processing).  The  conditional  Cramer  Rao  lower  bound  (CRLB)  is  developed  for  the  Direction-Of- 
Arrival  (DOA)  estimation  of  two  unresolved  Rayleigh  targets  using  a  standard  monopulse  radar. 
Then  the  modified  CRLB  is  used  to  give  insight  into  the  boresight  pointing  for  monopulse  DOA 
estimation  of  two  unresolved  targets.  Monopulse  processing  is  considered  for  DOA  estimation  of 
two  unresolved  Rayleigh  targets  with  known  or  estimated  relative  Radar  Cross  Section  (RCS). 
The  performance  of  the  DOA  estimator  is  studied  via  Monte  Carlo  simulations  and  compared 
to  the  modified  CRLB. 

Keywords:  Radar  Signal  Processing,  Target  Tracking,  Sensor  Resolution 
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1.  INTRODUCTION 


While  the  problem  of  tracking  multiple  targets  has  been  studied  extensively  in  recent  years, 
the  issue  of  finite  sensor  resolution  has  been  completely  ignored  in  almost  all  studies  [1].  Typi¬ 
cally,  the  targets  are  assumed  to  be  detected  with  a  given  probability  of  detection  in  the  presence 
of  false  alarms  and  clutter,  and  the  target  measurements  are  modeled  as  the  true  values  plus 
independent  Gaussian  errors  [2].  However,  when  two  targets  are  closely  spaced  with  regard  to 
the  resolution  of  the  sensor,  the  measurements  from  the  two  targets  can  be  merged  into  a  single 
measurement.  The  Joint  Probabilistic  Data  Association  (JPDA)  algorithm  [2]  was  extended 
in  [3]  to  develop  the  JPDAM  algorithm,  which  includes  possibly  merged  measurements  (i.e., 
unresolved  targets).  Multiple  Hypothesis  Tracking  (MHT)  [4]  was  extended  in  [5]  to  include 
possibly  unresolved  measurements.  In  the  JPDAM  algorithm  or  MHT,  merged  measurements 
are  modeled  as  a  “power”  weighted  centroid  of  the  two  targets.  While  the  application  of  the 
results  of  this  paper  can  be  used  to  enhance  the  JPDAM  and  MHT  algorithms,  the  focus  of  this 
paper  is  the  generation  of  two  Direction-Of- Arrival  (DOA)  estimates  for  tracking  two  unresolved 
targets  rather  than  the  use  of  a  centroid  for  the  two  targets. 

Monopulse  is  a  simultaneous  lobing  technique  for  determining  the  angular  location  of  a  source 
of  radiation  or  of  a  “target”  that  reflects  part  of  the  energy  incident  upon  it  [6].  In  an  amplitude 
comparison  monopulse  radar  system,  a  pulse  of  energy  is  transmitted  directly  at  the  predicted 
position  of  the  target,  and  the  target  echo  is  received  with  two  beams  that  are  squinted  relative 
to  the  predicted  position  of  the  target.  (Note  that  two  beams  are  required  for  each  angular 
coordinate.)  Traditionally,  the  DOA  of  a  target  is  estimated  with  the  in-phase  part  (i.e.,  the  real 
part)  of  the  monopulse  ratio,  which  is  formed  by  dividing  the  difference  of  the  two  received  signals 
by  their  sum.  When  target  echoes  interfere  (i.e.,  the  echoes  are  not  resolved  in  the  frequency 
or  time  domains),  the  DOA  estimate  indicated  by  the  in-phase  monopulse  ratio  can  wander  far 
beyond  the  angular  separation  of  the  targets  [6].  In  [7],  Daum  studied  the  angular  estimation 
accuracies  of  a  standard  monopulse  radar  system  for  two  unresolved  targets  and  showed  that  the 
estimation  accuracies  decline  significantly  as  the  variance  of  the  target  amplitude  fluctuations 
increases.  A  generalized  maximum  likelihood  ratio  test  is  developed  in  [8]  for  detecting  the 
presence  of  unresolved  Rayleigh  targets.  While  the  problem  of  DOA  estimation  of  unresolved 
targets  has  been  addressed  in  several  studies  involving  array  signal  processing  [9,10]  and  multiple- 
beam  monopulse  (i.e.,  more  than  two  beams  per  angular  coordinate)  [11],  the  work  of  Sherman 
in  [6,12]  is  the  only  technique  that  utilizes  a  standard  monopulse  system  in  the  DOA  estimation 
of  two  unresolved  targets. 

Since  the  hardware  of  many  existing  systems  provides  output  signals  for  standard  monopulse, 
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techniques  which  directly  use  this  output  can  be  implemented  in  timely  and  cost-efficient  manner 
without  the  development  of  new  hardware.  Sherman  proposed  in  [6,12]  the  use  of  complex 
monopulse  ratios  from  two  pulses  separated  sufficiently  in  time  so  that  the  relative  phase  of 
the  two  targets  changes,  but  sufficiently  close  in  time  so  that  the  amplitudes  of  the  two  targets 
remain  fixed.  Sherman  then  used  the  five  measured  quantities  (i.e.,  the  in-phase  and  quadrature 
components  of  the  two  complex  monopulse  ratios  and  the  ratio  of  the  measured  amplitudes 
of  the  sum  signal)  to  compute  the  DOAs  of  the  two  targets,  the  two  relative  phases,  and  the 
ratio  of  the  target  amplitudes.  In  [12],  the  in-phase  and  quadrature  monopulse  ratios  and 
the  Method  of  Moments  (MM)  were  used  to  develop  estimators  of  the  DOAs  for  tracking  two 
unresolved,  fixed- amplitude  targets  with  known  relative  Radar  Cross  Section  (RCS).  Berkowitz 
and  Sherman  showed  in  [13]  that  under  similar  assumptions  the  in-phase  monopulse  ratios  from 
the  two  angular  coordinates  and  the  ratio  of  the  measured  amplitudes  from  two  pulses  can  be 
used  to  estimate  the  centroid  of  the  two  targets  and  the  slope  of  the  line  connecting  them. 
However,  achieving  two  pulses  with  echoes  that  satisfy  the  requirements  of  Sherman’s  technique 
is  questionable.  Furthermore,  Sherman  utilized  a  deterministic  formulation  of  the  problem  to 
develop  his  approach.  Thus,  to  date,  the  results  of  Sherman  have  not  been  further  developed 
and  reported  in  the  literature,  and  no  algorithm  is  available  for  estimating  the  DOAs  of  two 
unresolved  Rayleigh  targets. 

In  this  paper,  the  results  on  monopulse  statistics  derived  in  [8]  are  used  to  develop  estimators 
that  generate  two  DOA  estimates  to  be  used  in  tracking  rather  than  relying  on  the  centroid 
measurement  conventionally  used  in  tracking.  The  use  of  two  DOA  estimates  will  be  critical  for 
the  tracking  of  two  aircraft  flying  in  formation  or  separating  ballistic  missiles,  where  the  flight 
paths  result  in  the  targets  being  unresolved  for  more  than  a  couple  of  measurements.  The  two 
DOA  estimates  are  computed  under  the  assumption  of  a  known  or  estimated  relative  RCS  of  the 
two  Rayleigh  targets,  which  can  be  obtained  through  various  methods.  For  example,  a  priori 
information  about  the  targets  such  as  type  and  flight  pattern  could  be  used  to  estimate  the 
relative  RCS.  If  RCS  estimates  of  the  individual  resolved  targets  were  obtained  while  tracking, 
they  can  also  be  used.  Beam  or  spatial  agility  could  be  exploited  to  estimate  the  relative  RCS. 
However,  for  this  paper,  the  relative  RCS  is  treated  as  a  known  or  estimated  quantity. 

Some  background  material  and  notation  are  given  in  Section  2,  while  the  probability  distri¬ 
bution  of  the  monopulse  measurements  for  two  unresolved  Rayleigh  targets  is  given  in  Section 
3.  In  Section  4,  the  conditional  Cramer  Rao  Lower  Bounds  (CRLBs)  are  developed  for  the  DOA 
estimation  of  two  unresolved  Rayleigh  targets  using  a  standard  monopulse  radar.  The  modified 
CRLBs  (MCRLB)  [14],  which  are  needed  to  remove  the  dependence  on  the  received  amplitudes, 
are  calculated  and  then  used  to  study  the  effects  of  antenna  boresight  pointing  on  the  DOA 
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estimation.  The  DOA  estimation  for  the  case  of  two  unresolved  targets  with  known  relative 
RCS  is  treated  in  Section  5,  where  the  mean  of  the  in-phase  inonopulse  ratio  and  the  variance 
of  the  in-phase  and  quadrature  monopulse  ratios  are  utilized  to  estimate  the  DOAs  of  the  two 
targets.  Simulation  results  that  illustrate  the  performance  of  the  DOA  estimators  and  compares 
it  to  the  MCRLB  are  also  given  in  Section  5.  Concluding  remarks  are  given  in  Section  6. 

2.  DEFINITIONS  AND  BACKGROUND 

The  in-phase  and  quadrature  parts  of  the  sum  and  difference  signals  for  two  independent, 
unresolved  Rayleigh  targets  can  be  written  as 


si  =  c*i  cos  (p i  +  a2  cos  <p2  +  nsi  (1) 

sq  =  ai  sin  <p i  +  a2  sin  <p2  +  tisq  (2) 

di  =  aiTji  cos  (fri  +  a27y2  cos  (p2  -I-  ndi  (3) 

( Iq  =  airji  sin  (pi  +  a2r/2  sin  <p2  +  ndQ  (4) 


where  a;  is  the  Rayleigh  distributed  amplitude  of  the  received  signal  from  target  i,  (pi  is  phase 
of  the  received  signal  from  target  i,  rji  is  the  unitless  (volt/volt)  DOA  for  target  i ,  and 

nSi  ~  N(0,  <t|)  usq  ~  N( 0,  cr|) 

ndi  ~  N(0,  aj)  ndQ  ~  ^(0,  cr^) 

The  notation  N(x,a 2)  denotes  a  Gaussian  distribution  with  mean  x  and  variance  a2.  For  this 
paper,  the  noise  terms  nsi,nsQ,ndi:  and  ndQ,  are  assumed  to  be  independent. 

Letting  A  and  ip  denote  the  measured  amplitude  and  phase  of  the  sum  signal  gives 

si  =  A  cos  ip  sq  =  A  sin  ip  (5) 


Since  op  and  a2  are  Rayleigh  distributed  and  (p\  and  cp2  are  uniformly  distributed  on  (— 7r,  7r], 
si  and  Sq  are  Gaussian  random  variables.  Applying  the  transformation  of  random  variables  in 
(5)  to  the  PDF  of  s/  and  sq  gives  the  PDF  of  the  observed  SNR  as 


3Zr  +  1 


+ 


>  0 


Revised  manuscript  submitted  to  IEEE  Trans,  on  Aerospace  and  Electronic  Systems,  August  10,  2000 


4 


where  !ft r  is  the  SNR  parameter  of  the  Rayleigh  signal  given  by 

(8) 

where  E[- }  denotes  the  expected  value.  Since  the  relative  RCS  of  the  targets  is  assumed  to  be 
known,  let  ^Rr2  =  ASR^,  A  >  0.  Then  -ft#  =  (1  +  A)3R/ji-  For  N  subpulses  at  distinct  frequencies 
(i.e.,  independent),  the  Maximum  Likelihood  (ML)  estimate  of  is  given  by 


where  5R0fe  denotes  the  observed  SNR  for  subpulse  k.  Then  is  an  unbiased,  efficient  estimator 
of  5R  with  variance  given  by  VAR[$ft^|5Rft]  =  iV-1(5jR  +  l)2. 

Denoting  s  =  si  +jsQ  and  d  =  di+jdQ,  the  in-phase  and  quadrature  parts  of  the  monopulse 
ratio  are  given  by 

,d x  sjdi  +  sQdQ 

VI  =  Reb>  =  -if  +4-  (10) 

Mdx  sidQ-SQdi 

!>=  +  4  (11) 


3.  STATISTICS  OF  THE  MONOPULSE  RATIOS 

The  PDF  of  the  monopulse  ratios,  yi  and  yQ,  for  a  single  pulse  can  be  obtained  by  application 
of  the  transformation  of  random  variables  in  (5),  (10)  and  (11)  as  in  [8]  to  /(s/,  sq,  d/,  dg|$), 
where  #  denotes  the  parameter  set  {5Jri,  3Jr2,  VhV2i  &d}-  Integrating  the  result  with  respect 
to  ip  and  conditioning  the  density  on  A  in  the  form  of  gives  [8] 

f(yi,  yQ\fto,  $)  =  f(yiWo,  *)f(yQ  |»o,  *)  (12) 

where 

/(wl»«.#)  =  jv(g,,^-) 

/(*,!*.,*) 

_  _  ^RlVl  +  ^R2V2 
yi  +  ^R2  +  1 

=  rcr2  Kmtf  +  XR2vt  +  XRiKR2(Vi-ri2)2 

■<rl  Kri  +  5Rr2  +  1 


(13) 

(14) 

(15) 

(16) 
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Thus,  yi  and  yq  are  conditionally  independent,  Gaussian  random  variables  with  a  common 
variance.  The  mean  of  yi  is  a  “power”  weighted  average  of  the  DOAs  of  the  two  targets  and 
independent  on  5R0,  while  the  mean  of  yQ  is  zero. 

Using  N  subpulses  at  distinct  frequencies  to  ensure  that  the  targets’  amplitudes  are  indepen¬ 
dent  for  each  subpulse,  the  conditional  ML  estimate  of  yj  given  !R0  for  each  pulse  is 


N  _i  N 

Vi  =  |  ^  ^  3?onJ  ^  '  ^okVlk 

n  =  l  k= 1 


1 

NY, n 


N 

^  ]  5? okyik 

k= 1 


(17) 


where  and  yik  denote  the  observed  SNR  and  in-phase  monopulse  ratio  for  subpulse  k  and 
Ypf  is  given  by  (9).  Note  that  while  the  ML  estimate  of  yj  is  conditional,  yi  is  independent  of 
by  (15).  Since  the  yik  are  Gaussian  random  variables,  yj  is  the  minimum  variance  estimate 
of  yi  and  a  Gaussian  random  variable  with  variance  given  by 


yi 


N 


-i  -1 


u=i 


Q 

2  NYn 


(18) 


Then  <7?  is  the  conditional  CRLB  for  any  unbiased  estimator  of  yi  given  { T?0fc } i  -  The 
term  conditional  here  is  used  to  denote  the  fact  that  the  Fisher  information  and  CRLB  are 
developed  with  the  amplitude-conditioned  PDF.  While  cr?;  provides  the  variance  of  yj  for  real¬ 
time  or  actual  tracking,  cr?7  cannot  be  used  for  performance  prediction  because  it  is  a  function 
of  {$ok}k=v  which  are  measured  quantities.  The  MCRLB  [14]  is  obtained  by  averaging  the 
conditional  Fisher  information  with  respect  to  the  measured  amplitudes  of  the  received  pulses. 
While  the  MCRLB  is  a  somewhat  looser  bound  than  the  CRLB,  it  can  be  used  for  performance 
prediction.  Since  /(yj|5ft0,  <f>)  satisfies  the  regularity  conditions  with  respect  to  yi  [14],  the 
MCRLB  is  a  lower  bound  on  the  variance  of  any  unbiased  estimator  of  yi ,  which  is  given  by 

-  Ewr' = 2mm  +%-Ti)  (19> 

where  yt  is  any  unbiased  estimator  of  y/. 

Figure  1  gives  a  comparison  of  the  MCRLBs  of  variance  of  yj  for  a  single  target  and  two 
unresolved,  equal-amplitude  targets  versus  the  total  SNR  in  a  single  frequency  (i.e.,  N  =  1).  The 
MCRLB  of  the  variance  of  yz  for  a  single  target  is  obtained  by  setting  5R^2  =  0  in  (19).  For  the 
two-target  case,  the  total  SNR  is  given  by  3?# l  +  In  order  to  obtain  yi  —  0  in  both 

cases,  the  single  target  was  set  at  the  boresight,  and  the  two  targets  were  situated  symmetrically 
about  the  boresight  (i.e.,  rj\  =  —772)  and  separated  by  one-half  of  a  beamwidth  with  r]\  =  0.4. 
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Figure  1.  MCRLB  of  the  variance  of  i/j  for  a  single  pulse  and  single  target  at  the  boresight  and 
two  targets  separated  by  one-half  beamwidth  and  symmetric  about  the  boresight 


Figure  1  shows  that  a  total  SNR  of  14  dB  gives  a  MCRLB  of  0.02  for  the  single-target  case  and 
0.095  for  the  two-target  case.  Figure  1  also  shows  that  doubling  the  energy  in  a  single-frequency 
waveform  to  obtain  a  total  SNR  of  17  dB  gives  a  MCRLB  of  0.01  for  the  single-target  case  and 
0.09  for  the  two-target  case.  Thus,  doubling  the  energy  in  a  single-frequency  pulse  gives  only  a 
small  reduction  in  the  MCRLB  for  the  two-target  case.  However,  if  the  energy  in  the  pulse  was 
doubled  by  adding  a  second  frequency  so  that  two  independent  observations  of  the  unresolved 
targets  is  obtained,  the  variance  would  be  reduced  by  50  percent  because  the  errors  are  Gaussian 
(not  shown  in  figure).  Therefore,  frequency  agility  is  critical  to  improving  the  monopulse  angle 
estimation  when  two  unresolved  targets  are  present. 

4.  CONDITIONAL  AND  MODIFIED  CRLBS  FOR  DOA  ESTIMATION 

The  conditional  and  modified  CRLBs  associated  with  r)\  and  r?2  are  developed  in  this  section. 
The  MCRLB  is  used  here  to  study  the  effects  of  antenna  pointing  on  the  DOA  estimation  and 
it  is  used  in  Section  5  to  assess  the  performance  of  the  DOA  estimators. 

The  conditional  Fisher  Information  Matrix  (FIM)  [15,  p.  79]  associated  with  fj\  and  772  based 
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on  N  observations  of  yi  and  yq  is  given  by 


Iy„yQ(m,V2\N,  0)  = 


(20) 


N 


-E 

k=i 


E 


E 


02lnf(yIk,yQk\^ok-^) 


dyj 

rg2Iii/(y/fc;yQfc|^ofc,<l>) 


dvidm 


$ok,$ 


E 

E 


d2in/(y/fc)  yQfcl^oA^  jO 
drjidr)2 

^2ln/(j//fci  yQfcjjftofc)  $) 

dy\ 


where  f(yikl,  yQkffiok,  $)  is  the  conditional  PDF  of  the  monopulse  ratios  for  pulse  k.  Using  (12) 
through  (15)  in  (20)  gives 


Iyi,yQ(f)ufl2\N,{Xok}Z=1,$) 

2 


2NtoR19tR2 


X 


Q{^R1  +  ^R2  +  l)2 

»jur~  2/  „  .  \2i  ..  2 


(21) 


Rr2 
2 

Q 


YjV  +  -  ^1  +  ]  Yn  +  -  +  $lR2&vj  (jl2  -  ^Rl^V^j 

SRr2 


Yn  +  -  (jli  +  ^R2Ayj  (r/2  -  5Ri?iAr/^  Yn  +  -  (y2  -  ytmAy'j 


where  Ay  =  771  —  y2  >  0.  The  dependence  of  Iy,,yQ{fh,  772 1  AT,  {K0fc}fcLi>4>)  on  Yn  can  removed 
by  averaging  over  the  observed  SNR.  {JRofc}^,  thereby  deriving  the  modified  FIM, 


(22) 


which  is  equivalent  to  setting  Yn  —  i?[y)v|$]  =  5Rfii  +  9?k2  +  1  in  (21)  since  Yn  appears  only 
linearly. 

A  bound  on  the  covariance  of  any  unbiased  estimator  of  771  and  772  [15,  p.  79]  is  given  by  the 
conditional  CRLB, 


COV 


m 

V2 


> 


Jyi,yQ(rh,fl2\N,{Mok}Z=l^)=  [lv,,yQm,  fi2\N,{nok}%=v*)]  1  (23) 


where 


Jyi,yQ{m^2\N,  {RofclfcU.S)  = 


X 


-1  - 


^R2 
$Rl 
2 


4iV5RR1!Rii2A7?2 

2n 


(24) 


qYN 


[I  +  iiv(,'2's',lA’')  1 

(r?i  +  5R/j2A?7^  (772  -  3?RiA77^ 


-1  - 


qYN 
3?R2 


(pi  +  $R.2&v)  (v2  ~  tomAqJ 


1  + 


qYn 


(vi  +  ^■R2^r/Sj 


In  (23),  the  matrix  inequality  for  matrices  A  and  B  with  A  >  B  is  defined  as  A  —  B  is  nonnegative 
definite.  Since  f(yi,  t/q|5R0,  <E>)  satisfies  the  regularity  conditions  with  respect  to  771  and  772  [14], 
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the  MCRLB  is  a  lower  bound  on  the  covariance  of  any  unbiased  estimator  of  rfo  and  772  [14]. 
The  MCRLB,  Jyi,yQ(fji,fh\N> $),  is  derived  by  using  the  modified  FIM  of  (22)  instead  of  the 
conditional  FIM  of  (21)  in  (23).  Thus,  the  MCRLB,  jy,,yQ{viiV2\N,<fr),  is  obtained  by  setting 
Yn  =  #[yjv|3>]  =  +  ?R.r2  +  1  in  (24)  and  given  by 


4N^Rri^.R2^V2 

,  2(772  -  ^iAT?)2 


$R2  f,  2(772  -  ytjuArjy 

I-  +  $?R2  +  1)- 

_  2(771  +  jftfl2A77)(772  -  UriAtj) 

q($Ri  +  +  1) 


_  2(r/l  +  ^R2^n)(j]2  -  ^RiAt?) 

qffim  +  ^R2  + 1) 

L  2(771  +  ^2^77) 2  I 
3^772  I-  q(^-Rl  +  ^R2  +  1)  -I 


Although  the  MCRLB  is  known  to  be  a  somewhat  looser  bound  than  the  CRLB,  it  is  needed 
to  remove  the  dependence  on  the  received  amplitudes  so  that  performance  prediction  can  be 
accomplished.  The  diagonal  elements  of  Jyi,yQ (t?i,  772I./V,  <F)  give  the  MCRLBs  for  the  variances 
of  7)1  and  7)2  and  are  used  to  assess  the  efficiency  [15,  p.  71]  of  the  DOA  estimators  in  the  next 
section. 

The  MCRLB  of  (25)  can  be  used  to  study  the  effects  of  sensor  pointing  on  the  DOA  estimation. 
For  this  study,  the  effects  of  the  antenna  gain  pattern  are  not  assumed  to  be  included  in  A,  the 
relative  RCS  of  the  two  targets.  The  effects  of  the  antenna  gain  pattern  can  instead  be  included 
in  the  anaylsis  of  the  MCRLB  by  using 

Sm  =  ^Ricos4(^-)  (26) 

^2  =  ^r2cos4(-^)  (27) 

in  (25)  for  5ft/*i  and  5?r2i  respectively.  The  r]bw  denotes  the  DOA  value  at  the  one-way,  half-power 
point  on  the  antenna  gain  pattern  of  the  sum  channel.  Thus,  at  rji  =  r]bw ,  3?.ri  =  —  6  dB. 

For  a  monopulse  error  slope  km  in  beamwidths,  2r]bw  ~  km.  For  all  of  the  examples  in  this  paper, 

T]bw  —0.8. 


Figure  2(a)  shows  the  average  MCRLB  for  two  19  dB  targets  (i.e.,  5ft/u  =  5ft .R2  =  19  dB  in  (26) 
and  (27)  for  N  =  1)  separated  by  one-half  beamwidth  (i.e.,  A77  =  rjbw  =  0.8)  versus  the  DOA 
of  target  1,  771,  for  N  =  1,  2,  4,  8  and  16  subpulses.  The  average  MCRLB  is  the  average  of  the 
two  diagonal  elements  of  (25).  For  each  case  in  Figure  2(a),  the  energy  in  the  transmitted  pulse 
is  fixed  and  it  is  divided  into  N  subpulses  at  different  frequencies.  Thus,  for  N  —  8  subpulses, 
5ft^x  =  5fti*2  =  10  dB  in  (26)  and  (27).  The  DOA  for  target  2  is  given  by  772  =  771  -  Arj. 
Thus,  771  =  0  in  Figure  2(a)  corresponds  to  target  1  on  the  antenna  boresight,  while  r]i  =  0.8 
corresponds  to  target  2  on  the  antenna  boresight.  Also,  r/i  —  —0.2  and  772  =  —1.0  in  Figure  2(a) 
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corresponds  to  pointing  the  antenna  boresight  on  the  outside  of  the  two  targets.  The  average 
of  the  two  MCRLBs  in  Figure  2(a)  shows  that  pointing  the  antenna  boresight  exactly  between 
the  two  targets  (i.e.,  77 \  =  0.4)  gives  the  smallest  average  MCRLB  for  N  =  4.  8,  and  16  and  the 
largest  bound  for  pointing  between  the  targets  when  N  =  1.  For  N  =  1,  the  minimum  average 
MCRLB  is  achieved  by  pointing  the  antenna  boresight  just  to  the  outside  of  one  of  the  two 
targets.  Figure  2(a)  also  shows  that  the  average  MCRLB  is  rather  insensitive  to  the  antenna 
boresight  when  it  is  positioned  between  the  two  targets.  For  N  =  16,  the  rapid  increases  in 
the  average  MCRLB  near  the  edges  of  the  plot  are  due  to  the  rapid  increase  of  the  MCRLB 
for  the  weaker  target.  To  illustrate  this  point,  Figure  2(b)  gives  the  individual  MCRLBs  of  the 
two  targets  for  N  =  4  and  8.  Figure  2(b)  shows  that  pointing  the  antenna  boresight  closer  to 
target  1  gives  a  smaller  MCRLB  for  the  variance  of  r)i  and  a  larger  MCRLB  for  the  variance 
of  7)2.  Note  the  rapid  increase  in  the  MCRLB  of  the  weaker  target  when  the  antenna  is  pointed 
to  the  outside  of  one  of  the  two  targets.  Also,  note  that  the  reduction  in  the  MCRLB  gained 
through  frequency  agility  (i.e.,  N  =  8  rather  than  N  =  4)  is  lost  due  to  the  low  subpulse  SNR 
at  r)i  =  -0.15  for  772  and  r/i  =  0.95  for  771. 

If  the  DOA  estimate  for  a  single  target  is  desired,  Figure  2(b)  shows  that  pointing  the  antenna 
closer  to  or  beyond  the  target  of  interest  (but  away  from  the  interfering  target)  minimizes  the 
MCRLB  for  that  target.  The  minimum  MCRLB  for  the  variance  of  7)1  with  N  —  4  is  0.015 
at  Tj\  =  —0.3  and  with  N  =  8  is  0.012  at  771  =  -0.2.  For  N  =  1  (not  shown  in  Figure  2(b)), 
the  minimum  MCLRB  for  the  variance  of  171  is  0.012  at  rji  =  —0.5.  For  the  single  target  DOA 
estimation,  an  optimal  pointing  angle  always  exists,  which  moves  towards  the  target  of  interest 
as  N  increases  for  the  given  separation  of  the  two  targets  and  target  amplitudes.  For  N  >  1, 
note  that  pointing  the  antenna  boresight  away  from  the  intefering  target  is  optimal  for  single 
DOA  estimation,  whereas  pointing  between  the  targets  is  optimal  when  estimating  both  DOAs. 

Figure  2(c)  gives  further  insight  into  the  pointing  of  the  antenna  boresight  for  the  DOA 
estimation  of  two  unresolved  targets.  Using  eight  frequencies  (i.e.,  IV  =  8)  to  simultaneously 
estimate  both  DOAs  gives  a  minimum  MCRLBs  for  the  individual  targets  of  0.028  with  the 
antenna  boresight  pointed  exactly  between  the  two  targets,  as  seen  also  in  Figure  2(a).  If  the 
eight  subpulses  are  used  in  two  dwells  of  four  subpulses  each,  the  two  DOAs  can  be  sequentially 
estimated  by  using  the  first  dwell  to  estimate  the  nearer  target  DOA,  and  using  the  second  to 
estimate  the  other  target  DOA.  This  method  also  gives  a  minimum  MCRLB  of  0.028  obtained 
by  pointing  to  the  outside  of  the  targets.  Further  analysis  showed  that  if  the  separation  between 
the  targets  is  larger  than  one  half  of  the  beamwidth  (i.e.,  A77  >  775^ ) ,  sequential  estimation  yields 
a  lower  MCRLB,  whereas  if  the  separation  is  less  than  one  half  the  beamwidth,  simultaneous 
estimation  provides  the  best  performance  bound  (not  shown).  Therefore,  for  equal  amplitude 


Revised  manuscript  submitted  to  IEEE  Trans,  on  Aerospace  and  Electronic  Systems,  August  10,  2000 


10 


targets,  the  case  A rj  =  t?{,w  shown  in  Figure  2(c)  denotes  a  boundary  in  the  decision  process  for 
the  optimal  pointing  strategy. 

Figure  3(a)  shows  the  average  MCRLB  for  two  19  dB  targets  (i.e.,  =  19  dB  in 

(26)  and  (27)  for  N  =  1)  separated  by  one-fourth  beamwidth  (i.e.,  A77  =  0.5r}t,w  =  0.4)  versus 
rji  for  N  =  1,  2,  4,  8,  and  16  subpulses.  Thus,  r?i  =  0  in  Figure  3(a)  corresponds  to  target  1 
on  the  antenna  boresight,  while  771  =  0.4  corresponds  to  target  2  on  the  antenna  boresight.  The 
average  of  the  two  MCRLBs  in  Figure  3(a)  shows  that  pointing  the  antenna  boresight  exactly 
between  the  two  targets  (i.e.,  771  =  0.2)  gives  the  smallest  average  MCRLB  for  all  values  of  AT. 
Figure  3(a)  also  shows  that  increasing  the  number  of  frequencies  from  N  =  8  to  N  =  16  actually 
results  in  an  increase  in  the  average  MCRLB.  Note  that  this  result  differs  from  that  shown  in 
Figure  2(a),  and  it  is  thought  to  be  the  result  of  a  low  subpulse  SNR  at  N  =  16  and  a  lower 
average  MCRLB  in  3(a)  compared  to  that  in  2(a). 

Figure  3(b)  gives  insight  into  the  pointing  of  the  antenna  boresight  for  the  DOA  estimation 
of  a  single  target  in  the  presence  of  another  target  for  one-fourth  beamwidth  separation.  It  is 
interesting  to  note  the  minimum  MCRLB  for  single  DOA  estimation  corresponds  to  pointing 
between  the  two  equal  amplitude  targets  for  N  =  8,  in  contrast  to  the  results  of  Figure  2(b). 
A  comparision  of  the  sequential  and  the  simultaneous  estimation  algorithms  for  this  case  in  a 
manner  similar  to  that  of  Figure  2(c)  indicates  that  it  is  always  advantageous  to  use  simultaneous 
estimation. 

Figure  4(a)  shows  the  average  MCRLB  for  a  19  dB  target  and  a  16  dB  target  (i.e.,  Ji/jj  = 
25Ri?2  =  19  dB)  separated  by  one-fourth  beamwidth  versus  rji  for  N  =  1,  2,  4  and  8  subpidses. 
For  the  N  —  1  and  2  cases,  the  optimal  pointing  strategy  is  to  steer  the  antenna  boresight 
slightly  towards  the  stronger  target,  while  for  the  other  cases,  steering  should  lean  toward  the 
weaker  target.  It  appears  that  the  simultaneous  DOA  estimator  favors  pointing  towards  the 
lower  power  target  when  subpulse  powers  are  too  low  ( N  =  4  or  8),  perhaps  to  avoid  very 
low  SNR  observations.  Note  also  that  the  advantages  of  frequency  agility  have  been  achieved 
by  N  =  4  and  the  degradation  in  performance  due  to  low  subpulse  SNR  can  be  seen  by  the 
steepening  of  the  graph  at  extreme  values  of  771  for  N  =  8.  The  MCRLB  for  the  single  target 
DOA  estimator  shown  in  Figure  4(b)  agrees  with  the  results  for  the  equal  energy  case  of  Figure 
3(b)  that  for  a  small  number  of  subpulses  the  antenna  should  be  pointed  at  or  outside  the 
position  of  the  targets  and  for  a  larger  number  of  subpulses  it  should  point  on  the  inside.  For 
the  case  of  two  unequal  amplitude  targets,  the  pointing  strategy  for  sequential  estimation  is  less 
clear.  Figure  4(b)  suggests  a  potential  iterative  approach.  The  pointing  of  the  first  dwell  is  in 
the  neighborhood  of  the  stronger  target  in  order  to  minimize  the  interference  from  the  weaker 
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Figure  2(a)  Average  MCRLB 


for  the  variances  of  fji  and  172  with  Ar;  =  77^  and  two  19  dB 


targets 


DIRECTION-OF-ARRIVAL  (DOA)  OF  TARGET  1 


Figure  2(b)  MCRLBs  for  variances  of  7)1  and  7)2  with  A77  =  r]^  and  two  19  dB  targets 
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Figure  2(c)  MCRLBs  for  sequential  and  simultaneous  estimation  of  r]i  and  772  with  A77  =  rjbw, 
two  19  dB  targets,  and  eight  subpulses  of  equal  energy  (Sequential:  N  =  4  in  two  consecutive 
dwells;  simultaneous:  N  —  8  in  one  dwell) 


target.  The  pointing  of  the  second  dwell  is  then  chosen  to  minimize  the  bound  in  estimating 
the  DOA  of  the  second,  weaker,  target  using  the  information  in  the  first  dwell. 

5.  DOA  ESTIMATION  FOR  TWO  UNRESOLVED  TARGETS 

As  seen  in  the  previous  section,  the  pointing  of  the  antenna  boresight  affects  the  information 
obtained  during  a  radar  dwell.  Based  on  the  results  of  Section  4,  two  sensor  pointing  strategies 
are  set  forth.  The  simultaneous  pointing  strategy  that  was  found  to  be  optimal  for  targets 
separated  by  less  than  one-half  beamwidth  involves  pointing  between  the  targets  and  estimating 
both  DOAs  from  data  gathered  in  one  dwell.  This  case  is  referred  to  as  the  fully  unresolved  case. 
The  second  strategy  proposed,  applicable  to  targets  which  are  separated  by  more  than  one-half 
of  a  beamwidth,  is  the  sequential  strategy  in  which  the  DOAs  are  estimated  from  data  gathered 
from  two  consecutive  dwells,  each  pointing  in  the  general  direction  of  one  target.  This  case 
can  be  thought  of  as  a  partially  resolved  target  scenario  since  estimation  of  each  parameter  can 
be  performed  individually.  In  this  section,  a  DOA  estimation  algorithm  based  on  the  Method 
of  Moments  (MM)  for  the  fully  unresolved  two-target  case  is  proposed.  The  partially  resolved 
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DIRECTION-OF-ARRIVAL  (DOA)  OF  TARGET  1 

Figure  3(a)  Average  MCRLB  for  the  variances  of  fji  and  %  with  At/  =  0.5 rjbw  and  two  19  dB 
targets 


Figure  3(b)  MCRLBs  for  variances  of  171  and  172  with  A77  =  0.5rjbw  and  two  19  dB  targets 
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Figure  4(a)  Average  MCRLB  for  the  variances  of  fft  and  7)2  with  Arj  =  0.5 f]bw  and  19  dB  and 
16  dB  targets,  respectively 


Figure  4(b)  MCRLB  for  the  variances  of  7)1  and  t/2  with  Arj  =  0.5%^  and  19  dB  and  16  dB 
targets,  respectively 
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target  case  is  briefly  discussed  in  [17]. 


For  the  case  of  two  fully  unresolved  Rayleigh  targets,  the  antenna  boresight  is  pointed  between 
the  targets  and  the  DOA  estimation  of  both  targets  is  accomplished  simultaneously  via  MM. 
For  the  MM,  two  sample  moments  are  required  for  the  estimation  of  two  DOAs.  Since  the 
observations  yik  are  conditionally  nonstationary  (i.e.,  the  variance  of  yjk  given  5?0/t  changes 
with  5R0fc  -  see  (13))  sample  moments  cannot  be  used  directly.  Thus,  yj  from  (17)  will  be  used 
to  obtain  the  first  expression  as 


^RlVl  +  ^R2V2 

IJ  T  ^  - 

Rftl  +  ^R2  +  1 


=  til 


(28) 


Let  Ar?  =  r/i  -  7j2  >  0,  so  that  r/i  >  q2-  Then,  for  »  1,  since  3?h2  =  AK/u,  where  A  >  0, 


A  A 

m  ~ yi  +  y+x^71 

1  A 

m  ~  yi  - 


(29) 

(30) 


Since  the  t/jfc  are  conditionally  nonstationary,  Gaussian  random  variables,  the  second  expres¬ 
sion  for  DOA  estimation  will  be  obtained  by  forming  a  Chi-squared  random  variable  with  2N-1 
degrees  of  freedom  from  yik,  K0fc,  and  yi  and  setting  the  random  variable  equal  to  its  mean. 
Thus,  for  N  >  1,  let 

Xn  =  [  Vl\  -  til  •••  VlN  -  til  VQ\  •••  UQN ]  (31) 

Rjv  =  2  diag[5i0i  ...  5R0tv  ^oi  •••  5R0/v]  (32) 

Since  for  N  =  1,  yi  =  yn,  let  X\  —  yQ\  and  R\  =  25R0i-  Also  let 

vN  =  XjjRNXNq-1  (33) 

where  q  is  given  by  (16).  Since  the  conditional  PDF  of  the  in-phase  and  quadrature  monopulse 
ratios  are  Gaussian  with  variance  given  by  (13)  through  (16),  v jv  is  a  Chi-squared  random 
variable  with  2 N  -  1  degrees  of  freedom  as  one  degree  of  freedom  has  been  lost  due  to  the 
estimation  of  yi  [16].  Then 

E[vN\$,{Mok}%=1]  =  2N  -1  (34) 

VARM$,  {?U-}£ Li]  =  2(2 N  -  1)  (35) 


Setting  vn  equal  to  its  mean  gives 


2N  —  1 


N  >  1 


(36) 
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Since  for  Ar/  >1.0  the  targets  would  be  processed  as  partially  unresolved  targets,  assume  that 
Ar/  <  1.0.  Thus,  pointing  the  antenna  boresight  between  the  targets  gives  |r/i|  <  1  and  |r/2|  <  1, 
which  allows  the  approximation  of  q  by 


.  9 

-T  +  TTTA7? 
(Tc  1  +  A 


(37) 


because  5Rri7/2  +  5ft_R2ftl  <<  ^Ri^R2{vi  -»/ 2)2  in  (16).  Note  that  (37)  implies  that  any  estimate 
of  q  must  be  greater  than  cr^/ag  and  less  than  cr^/ag  +  A(1  +  A)-15ftRimax(A7/2).  Inserting  q 
from  (36)  for  q,  using  (37)  in  (29)  and  (30),  and  applying  the  restrictions  on  an  estimate  of  q 
gives  DOA  estimates  as 


Vi  =  yi  + 


(38) 


V2 


(39) 


where 


(°, 

(1  +  A)5 


&r,  q  > 


q<°das2 

^Vbw^_ 
t  (1  +  A)2 
otherwise 


2-2 


Mr  +  Gd°s 


(40) 


t %w  is  defined  in  conjunction  with  (26)  and  (27),  5ft#  and  5ft  r  are  given  by  (8)  and  (9).  The  first 
case  of  (40)  is  introduced  to  ensure  that  the  DOA  estimates  are  real  numbers,  while  the  second 
is  introduced  to  prevent  the  difference  of  the  two  DOA  estimates  from  being  larger  than  2r)bw. 
Limiting  the  difference  of  the  two  DOA  estimates  was  found  to  be  critical  when  correcting  the 
known  relative  RCS  A  for  the  effects  of  the  antenna  gain  pattern  of  the  sum  channel. 


In  order  to  use  the  DOA  estimates  for  tracking,  estimates  of  the  variances  of  r)i  and  7)2  are 
needed.  As  detailed  in  the  Appendix,  the  variances  of  the  DOA  estimates  for  q  >>  cr^cr g2  can 
be  approximated  as 


VAR[t/2|<M  »ofc}*Li]*0 


1  A  q 

2iVyjv  +  25ft*(2 N  -  1  ){q-  afcf2)  _ 

-J—  +  _. _ £ _ 

2 NYn  2\Ur(2N  -  1  ){q-  o\a-g2) 


(41) 

(42) 


Setting  q  =  q  in  (41)  and  (42)  and  by  analysis  of  simulation  results  q(q  -  a^as2)  1  =  3.0  was 
found  to  provide  relatively  good  estimates  of  the  variances  for  the  DOA  estimates,  which  are 
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given  by 


&l  =  <? 


1 


2NYn 

1 


+ 


1.5A 


+ 


Kr(2N-1) 
1.5 


2NYn  YRr(2 N  -  1) 


(43) 

(44) 


Note  that  these  estimates  of  the  variances  are  functions  of  the  measured  data  and  cannot  be 
used  for  performance  prediction  of  the  estimators. 


Monte  Carlo  simulations  with  40000  experiments  were  conducted  to  study  the  performance 
of  the  DOA  estimators  for  various  values  of  N,  A 77,  and  A.  While  the  relative  RCS  of  the  two 
targets  is  assumed  to  be  known,  the  effect  of  the  antenna  gain  pattern  is  not  assumed  to  be 
included  in  A.  The  effects  of  the  antenna  gain  pattern  were  included  in  the  simulation  of  (l)-(4) 
by  using 


/  77i  7T  \  ,  , 

<*1  =  aqcos  - - )  (45) 

V4  rjbwJ 

q.2  =  a2cosf-^^-)  (46) 

V  4  T)bw  > 

for  a  1  and  02,  respectively.  The  effects  of  the  antenna  gain  pattern  were  addressed  in  the  DOA 
estimation  by  using  a  modified  A  in  (38)  through  (40),  which  is  given  by 


A  = 


(47) 


where  rj\  and  fj 2  are  the  DOA  estimates  that  result  from  ignoring  the  effects  of  the  antenna  gain 
pattern.  The  multiplier  modification  to  A  was  restricted  to  greater  than  0.25  and  less  than  4. 
Also,  during  the  simulations,  only  measurements  with  NY^  >  20  dB  were  utilized  in  the  DOA 
estimation  and  r]bw  =  0.8,  the  same  value  as  that  used  in  Section  4.  Analysis  of  the  simulation 
results  showed  that  the  use  of  measurements  with  a  total  SNR  less  than  20  dB  can  cause  the 
estimation  errors  to  deviate  significantly  from  the  MCRLB  and  the  estimator  to  be  unreliable. 
This  requirement  for  a  rather  large  SNR  is  most  likely  due  to  the  lack  of  tracking  information 
such  as  predicted  antenna  gains  for  the  targets.  While  an  SNR  of  20  dB  is  a  rather  strict 
requirement,  previous  methods  such  as  super  resolution  may  require  as  much  as  40  dB  [18]. 


The  effect  of  N  on  the  DOA  estimation  was  studied  using  A77  =  0.4  and  setting  N$tR  1  = 
N$lR 2  =  22  dB  (i.e.,  without  the  effects  of  the  antenna  pattern).  The  sample  averages  of  the 
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Figure  5(a)  Averages  of  Errors  in  DOA  Estimates  of  Target  1  for  At/  =  0.5rjbW  and  A  =  1 


Figure  5(b)  Modified  CRLB  of  Target  1  for  A?/  =  0.5t]i,w  and  A  =  1 
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errors  in  the  DOA  estimates  for  target  1  are  shown  versus  the  DOA  of  target  1  in  Figure  5(a). 
The  DOA  for  target  2  is  given  by  r/2  =  // 1  —  A r/.  Thus,  r/i  =  0  in  Figure  5(a)  corresponds 
to  target  1  on  the  antenna  boresight,  while  rji  =  0.4  corresponds  to  target  2  on  the  antenna 
boresight.  Figure  5(a)  shows  that  the  sample  average  errors  are  rather  small  for  each  case  of  N 
and  have  the  smaller  magnitudes  when  the  antenna  boresight  is  pointed  between  the  two  targets 
(i.e.,  r]i  «  0.2).  The  sample  standard  deviations  of  the  errors  in  the  DOA  estimates  for  target  1 
normalized  by  the  square  root  of  the  MCRLB  are  shown  versus  the  DOA  of  target  1  in  Figure 
5(b)  by  the  solid  lines.  Figure  5(b)  shows  that  the  efficiency  of  the  MM  estimator  improves 
(i.e.,  sample  variance  approaches  the  MCRLB)  as  N  increases  from  4  to  12  and  is  best  when  the 
antenna  boresight  is  pointed  at  target  1  (i.e.,  r/i  =  0).  However,  the  MM  estimators  are  shown  to 
be  rather  inefficient  in  Figure  5(b).  In  order  to  assess  the  effects  of  the  estimation  of  the  antenna 
gain  pattern  effects  on  the  relative  RCS  on  the  efficiency  of  the  MM  estimators,  the  simulation 
studies  were  conducted  with  A  as  a  known  quantity  rather  than  an  estimated  quantity  (i.e.,  the 
true  values  of  rft  and  772  were  used  in  (47)).  For  A  known,  the  sample  standard  deviations  of  the 
errors  in  the  DOA  estimates  for  target  1  normalized  by  the  MCRLB  are  shown  versus  the  DOA 
of  target  1  in  Figure  5(b)  by  the  dashed  lines.  Figure  5(b)  shows  that  the  efficiency  of  the  MM 
estimators  are  degraded  rather  significantly  by  the  estimation  of  A.  Thus,  if  the  two  targets  are 
under  track,  the  target  state  estimates  should  be  used  to  predict  the  effects  of  the  antenna  gain 
pattern.  Figure  5(b)  in  conjunction  with  Figure  3(a)  also  shows  that  increasing  the  number  of 
subpulses  at  distinct  frequencies  in  a  radar  pulse  improves  the  accuracy  and  efficiency  of  the 
DOA  estimation. 

The  effect  of  A r/  on  DOA  estimation  was  studied  using  N  =  8  and  setting  =  5R.R2  =  13 
dB.  The  averages  and  standard  deviations  of  the  errors  in  the  DOA  estimates  for  target  1  are 
shown  in  Figures  6(a)  and  6(b)  for  various  DOAs  of  target  1  and  A77  =  0.2,  0.4,  0.6,  and  0.8.  The 
DOA  estimation  for  A77  =  0.8  is  significantly  degraded  w'hen  compared  to  that  for  A77  =  0.2 
or  0.4.  Thus,  when  two  targets  are  separated  by  about  one-half  of  the  one-way  beamwidth, 
sequential  DOA  estimation  with  two  consecutive  dwells  at  the  individual  targets  will  be  better 
than  the  simultaneous  DOA  estimation  with  a  single  dwell.  This  observation  agrees  with  those 
of  Section  4. 

The  effect  of  A  on  DOA  estimation  was  studied  using  N  =  8,  A77  =  0.4,  and  =  16  dB. 
The  sample  averages  and  standard  deviations  of  the  errors  in  the  DOA  estimates  for  targets  1 
and  2  arc  shown  in  Figures  7(a)  and  7(b)  for  the  positive  DOAs  of  target  1,  the  corresponding 
negative  DOAs  of  target  2,  and  A  =  1,  0.5,  and  0.25  (i.e.,  5R/?2  —  16,  13,  and  10  dB).  The  771  =  0 
in  Figures  7(a)  and  7(b)  corresponds  to  target  1  on  the  antenna  boresight  and  772  =  -0.4,  while 
772  —  0  corresponds  to  target  2  on  the  antenna  boresight  and  771  =  0.4.  The  accuracy  of  the  DOA 
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SAMPLE  STANDARD  DEVIATION  OF  ERROR  w  SAMPLE  AVERAGE  ERROR 


Averages  of  Errors  in  DOA  Estimates  of  Target  1  for  N  =  8  and  A  =  1 


Figure  6(b)  Sample  Standard  Deviations  of  Errors  in  DOA  Estimates  of  Target  1  for  N 
and  A  =  1 
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Figure  7(a)  Averages  of  Errors  in  DOA  Estimates  for  TV  —  8  and  Arj  =  O.St^ 


Figure  7(b)  Standard  Deviations  of  Errors  in  DOA  Estimates  for  N  —  8  and  A r/  =  0.5 rj^ 
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◄ - Target  2 
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SOLID:  LAMBDA=1.0 


Figure  7(c)  Sample  Standard  Deviations  of  Errors  in  DOA  Estimates  Normalized  by  the  Esti¬ 
mates  of  Standard  Deviations  for  TV  =  8  and  A 77  =  0.5i%w 


estimation  for  target  1  improves  as  decreases,  while  accuracy  of  the  DOA  estimation  for 
target  2  degrades  as  5R«2  decreases.  The  errors  were  also  normalized  by  the  standard  deviation 
estimates  of  (43)  and  (44)  and  the  sample  standard  deviations  of  these  errors  are  given  in  Figure 
7(c),  which  indicates  that  the  variances  estimates  of  (43)  and  (44)  are  reasonably  good.  Figure 
7(c)  also  shows  that  the  variance  estimates  are  best  for  a  target  at  the  antenna  boresight. 

To  assess  the  impacts  on  the  DOA  estimation  of  uncertainty  in  A,  Monte  Carlo  experiments 
were  conducted  using  incorrect  values  for  the  relative  RCS,  A.  The  percent  error  in  A  is  define 
as  (A  —  A) A-1.  The  simulations  were  conducted  for  two  16  dB  targets  separated  by  one-fourth 
beamwidth  and  situated  symmetrically  about  the  antenna  boresight.  The  results  are  shown  in 
Figure  8.  For  this  case  of  equal  RCSs  (i.e.,  A  =  1.0),  negative  values  of  the  percent  error  indicate 
that  5ft  i  is  thought  to  be  greater  than  9?2,  while  positive  values  of  the  percent  error  indicate  that 
5ft  2  is  thought  to  be  greater  than  .  The  average  errors  in  Figure  8(a)  show  that  the  bias  in  the 
DOA  estimate  is  of  the  same  sign  as  the  percent  error.  The  biases  in  the  DOA  estimates  are 
due  to  the  DOA  estimator  shifting  the  estimates  toward  the  target  that  is  thought  to  be  weaker, 
which  is  the  result  of  matching  the  data  to  the  assumed  parameters.  For  less  than  25%  error  in 
A,  the  biases  in  the  DOA  estimates  is  less  than  10%  for  this  case.  Figure  8(b)  shows  that  the 
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effects  of  errors  in  A  on  the  standard  deviations  of  the  DOA  estimates  are  very  small.  Note  that 
the  standard  deviations  of  the  errors  are  smaller  for  the  target  that  is  thought  to  be  stronger. 
The  results  of  Figure  8  suggest  that  the  DOA  estimator  is  rather  robust  to  uncertainty  in  A. 

6.  CONCLUDING  REMARKS 

The  conditional  CRLB  and  modified  CRLB  were  developed  for  monopulse  DOA  estimation 
of  two  unresolved  Rayleigh  targets.  The  MCRLB  was  used  to  investigate  the  effects  of  frequency 
diversity  and  antenna  borcsight  steering  on  the  DOA  estimation.  For  a  given  waveform  energy, 
the  accuracy  of  the  DOA  estimation  was  found  to  improve  as  the  number  of  distinct  frequencies 
increases  and  as  long  as  the  expected  subpulse  energy  remains  above  a  few  dB.  The  study 
revealed  two  sensor  pointing  strategies  for  the  DOA  estimation:  simultaneous  and  sequential. 
The  simultaneous  pointing  strategy  involves  pointing  between  the  targets  and  estimating  both 
DOAs  from  data  gathered  in  one  dwell,  while  the  sequential  pointing  strategy  involves  estimating 
the  DOAs  with  data  gathered  from  two  consecutive  dwells,  each  pointing  in  the  general  direction 
of  one  target.  While  the  relative  performance  of  the  two  pointing  stategies  depend  on  the  SNRs 
of  the  two  targets  and  the  number  of  subpulses,  the  simultaneous  pointing  strategy  was  found  to 
be  better  for  targets  separated  by  less  than  one-half  beamwidth,  while  the  sequential  pointing 
strategy  was  found  to  be  better  for  targets  that  are  separated  by  more  than  one-half  of  a 
beamwidth. 

An  estimation  algorithm  based  on  the  Method  of  Moments  (MM)  was  developed  for  the 
simultaneous  DOA  estimation  of  two  unresolved  Rayleigh  targets  with  known  relative  amplitude. 
In  the  DOA  estimation,  the  effects  of  the  antenna  gain  pattern  were  treated  as  unknown  and 
addressed  in  the  DOA  estimation  with  a  one-step  iteration  on  the  DOA  estimates.  In  other 
words,  the  initial  DOA  estimates  obtained  by  ignoring  the  effects  of  the  antenna  gain  were 
used  to  predict  the  effect  of  the  antenna  gain  on  the  relative  amplitude  of  the  two  targets 
for  DOA  estimation.  The  performances  of  the  DOA  estimators  were  studied  via  Monte  Carlo 
simulations  and  compared  to  the  MCRLB.  The  results  of  the  Monte  Carlo  simulations  confirmed 
the  predicted  benefits  of  frequency  diversity  and  showed  that  frequency  diversity  also  improves 
the  efficiency  [15]  of  the  estimators.  The  results  support  the  use  of  simultaneous  DOA  estimation 
when  the  targets  are  separated  by  less  than  one-half  of  a  beamwidth  and  sequential  DOA 
estimation  when  the  targets  are  separated  by  more  than  one-half  of  a  beamwidth. 

The  algorithms  developed  in  this  paper  have  two  practical  limitations:  the  uncertainty  in  the 
antenna  gains  at  the  two  DOAs  and  the  uncertainty  in  the  relative  RCS.  The  results  show  that 
the  efficiency  of  the  MM  estimators  is  degraded  rather  significantly  by  the  estimation  of  the 
effects  of  the  antenna  gain.  If  the  two  targets  are  under  track  prior  the  merged  measurement, 
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Figure  8(a)  Averages  of  Errors  in  DOA  Estimates  for  N  =  8  and  A77  =  0.5 r]bw  with  Error  in  A 


Figure  8(b)  Standard  Deviations  of  Errors  in  DOA  Estimates  for  N  =  8  and  A q  =  0.5'rjbw  with 
Error  in  A 
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the  target  state  estimates  can  be  used  to  predict  the  individual  DOAs  so  that  the  degradation 
due  to  the  DOA  dependent  antenna  gain  can  be  reduced.  Also,  the  tracks  of  the  two  targets 
can  be  used  to  estimate  the  relative  RCS  with  sufficient  accuracy.  Thus,  this  paper  provides 
practical  algorithms  for  DOA  estimation  from  unresolved  measurements  using  estimated  values 
for  the  required  parameters  (relative  RCS  and  antenna  gains)  that  can  be  obtained  from  the 
tracks  on  the  two  targets. 


Revised  manuscript  submitted  to  IEEE  Trans,  on  Aerospace  and  Electronic  Systems,  August  10,  2000 


26 


REFERENCES 


1.  F.  Daum,  “A  System  Approach  to  Multiple  Target  Tracking,”  in  Multitarget-Multisensor 
Tracking:  Applications  and  Advances.  Vol.  II,  edited  by  Y.  Bar-Shalom,  Artech  House,  Inc., 
Norwood,  MA,  1992. 

2.  Y.  Bar-Shalom  and  X.R.  Li,  Multitarget-Multisensor  Tracking:  Principles  and  Techniques, 
YBS  Publishers.  Box  U-157,  Storrs,  CT,  06296-3157,  1995. 

3.  K.C.  Chang  and  Y.  Bar-Shalom,  “Joint  Probabilistic  Data  Association  for  Multitarget 
Tracking  with  Possibly  Unresolved  Measurements  and  Maneuvers,”  IEEE  Trans,  on  Auto. 
Cont.,  July  1984,  pp.  585-594. 

4.  Blackman,  S.S.,  Multiple-Target  Tracking  with  Radar  Applications,  Artech  House,  Inc., 
Norwood,  MA,  1986. 

5.  W.  Koch  and  G.  Van  Keuk,  “Multiple  Hypothesis  Track  Maintenance  with  Possibly  Unre¬ 
solved  Measurements,”  IEEE  Trans.  Aero.  Elect.  Sys.,  July  1997,  pp.  883-892. 

6.  S.M.  Sherman.  Monopulse  Principles  and  Techniques,  Artech  House,  Inc.,  Dedham,  MA,  1984. 

7.  F.  E.  Daum,  “Angular  Estimation  Accuracy  for  Unresolved  Targets,”  Proc.  of  1987  American 
Control  Conference,  Minneapolis,  MN,  June  1987,  pp.  1135-36. 

8.  W.D.  Blair  and  M.  Brandt-Pearce,  “Unresolved  Rayleigh  Target  Detection  Using  Monopulse 
Measurements,”  IEEE  Trans.  Aero.  Elect.  Sys.,  April  1998,  pp.  543-552. 

9.  A.  A.  Ksienski  and  R.  B.  McGhee,  “A  Decision  Theoretic  Approach  to  the  Angular  Resolution 
and  Parameter  Estimation  Problem  for  Multiple  Targets,”  IEEE  Trans.  Aero.  Elect.  Sys.,  May 
1968,  pp.  443-455. 

10.  S.  Haykin,  J.  Litva,  and  T.J.  Shepherd,  Eds.,  Radar  Array  Processing,  Springer-Verlag,  New 
York,  NY,  1993. 

11.  P.  Z.  Peebles  and  R.  S.  Berkowitz,  “Multiple-Target  Monopulse  Radar  Processing  Tech¬ 
niques,”  IEEE  Trans.  Aero.  Elect.  Sys.,  Nov.  1968,  pp.  845-854. 

12.  S.  M.  Sherman,  “Complex  Indicated  Angles  Applied  To  Unresolved  Targets  and  Multipath,” 
IEEE  Trans.  Aero,  and  Elect.  Sys.,  Jan.  1971,  pp.  160-170. 

13.  R.  S.  Berkowitz  and  S.  M.  Sherman,  “Information  Derivable  from  Monopulse  Radar  Mea¬ 
surements  of  Two  Unresolved  Targets,”  IEEE  Trans.  Aero,  and  Elect.  Sys.,  Sept.  1971,  pp. 
1011-1013. 

14.  F.  Gini,  R.  Reggiannini,  and  U.  Mengali,  “The  Modified  Cramer-Rao  Bound  in  Vector 
Parameter  Estimation,”  IEEE  Trans,  on  Communications,  vol.  46,  No.  1,  pp.  52-60,  January 
1998. 

15.  H.  L.,  Van  Trees,  Detection,  Estimation,  and  Modulation  Theory:  Part  I,  John  Wiley  &;  Sons, 
Inc.,  1968. 

16.  Y.  Bar-Shalom  and  X.R.  Li,  Estimation  and  Tracking:  Principles,  Techniques,  and  Software, 
Artech  House,  Inc.,  Dedham,  MA,  1993.  (Reprinted  by  YBS  Publishering,  Box  U-157,  Storrs, 
CT,  06296-3157,  1998). 

17.  W.  D.  Blair,  G.A.  Watson,  and  M.  Brandt-Pearce,  “Monopulse  Tracking  of  Two  Unresolved 
Rayleigh  Targets,”  Proceedings  of  the  Signal  and  Data  Processing  of  Small  Targets  1997, 
SPIE’s  International  Symposium,  San  Diego,  CA,  July  1997. 


Revised  manuscript  submitted  to  IEEE  Trans,  on  Aerospace  and  Electronic  Systems,  August  10,  2000 


27 


18.  D.  K.  Barton,  Ed.,  Radar  Resolution  and  Multipath  Effects.  Artech  House,  Inc.,  Dedham, 
MA,  1975. 


Revised  manuscript  submitted  to  IEEE  Trans,  on  Aerospace  and  Electronic  Systems,  August  10,  2000 


28 


Appendix 

Derivation  of  DOA  Estimator  Variance 


This  appendix  details  the  derivation  of  the  estimate  of  the  variance  of  rji  given  in  (41).  The 
estimate  of  the  variance  of  7)2  given  in  (42)  is  derived  similarly.  Using  r)i  as  defined  in  (38)  and 
the  orthogonality  principle  [15]  gives 


VAR[thl«, {»<*}"-i]  =  VARfo,|»,  {«<*}£.,]  +  i-VAR[V9|4,  {»*}£., 

R 

A 


+ 


2  NYn  ^ 


£[«]  -  I®.  {SMLiJ  +  {»0»}£..]](a.i) 


where  the  first  term  is  given  by  (18)  and  4qlwXRR{\  +  A)  2  >  q  —  cr^cr^2  >  0  from  (40).  From 
(33)  through  (36), 


Em,  {».»}£■]  =  9 

VARBI*.  {su}Ld  = 


(A.  2) 
(-4-3) 


Note  that  for  a  twice  differentiable  function  g(x), 


E[g{x)j  «  g(x )  +  1 g"(x)o2x  (A. 4) 

where  x  =  E[x],  =  VAR[x],  and  g"{x )  is  the  second  derivative  of  g(x)  with  respect  to  x 

evaluated  at  x  =  x.  Using  (A. 2)  through  (A. 4)  gives 

E[y/q-<rj<Ts2\^  OWfcU]  ~  \fq~  cr^S  2 - T  ( A .5) 

4{2N-l)yJ{q-ala-s2f 

Using  (A. 2)  and  (A. 5)  in  (A.l)  gives 


VAR[»H|*,{*oJfc}*U] 


l  _ V _ 

2  NYn  2&r(2N  -  1  )(q-  ajaf) 


8(2N  —  l)(q-  a2das2)i 


(A.6) 


For  q  —  cr^as2  >  0.5,  which  is  satisfied  by  any  nondegenerate  case  (i.e.,  rfo  772)  of  two  unresolved 
targets,  the  last  term  of  (A. 6)  can  be  ignored  to  obtain  (41). 
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ABSTRACT 

A  new  method  of  track  management  for  a  phased  array  radar  is  proposed  to  simplify  the  data  association  and 
improve  the  allocation  of  radar  resources.  The  tracks  are  organized  into  Association  Groups  and  Dwell  Groups. 
Association  Groups  are  used  for  course  gating  when  associating  new  measurement  data  with  existing  tracks,  while 
Dwell  Groups  are  used  to  efficiently  schedule  the  next  dwell  on  closely-spaced  targets.  By  considering  only  closely- 
spaced  tracks  in  the  same  Association  Group,  a  form  of  coarse  gating  is  inherently  done  to  cull  candidate  tracks  that 
are  unlikely  to  associate  with  a  measurement  from  a  dwell  directed  by  one  of  the  members  of  the  group.  A  Dwell 
Group  contains  tracks  that  are  spaced  sufficiently  to  allow  one  dwell  to  illuminate  all  of  its  members.  Dwell  groups 
lay  the  foundation  for  more  systematic  approaches  to  optimal  allocation  of  the  radar  resources.  Simulation  results 
are  presented  to  illustrate  the  new  technique. 

Keywords:  Multitarget  Tracking,  Radar  Signal  Processing,  Phased  Array,  Data  Association 

1.  INTRODUCTION 

In  the  past  radar  trackers  have  been  designed  with  a  keen  eye  on  the  computational  throughput  of  the  computers 
available.  The  designers  justifiably  traded  radar  resources  for  computational  simplicity,1  ,2  ,3  ,4  and.5  However, 
advances  in  computational  hardware  now  make  more  sophisticated  tracking  algorithms  feasible.  A  new  method  of 
track  management  for  a  phased  array  radar  is  proposed  to  make  more  efficient  use  of  radar  resources  at  the  expense 
of  added  algorithm  complexity  and  computational  load. 

The  approach  is  based  on  classifying  tracks  into  categories  depending  on  one  of  two  radar  tasks-  processing 
collected  data  or  gathering  new  data.  When  new  data  has  been  collected,  the  measurements  must  be  assigned  to 
tracks  according  to  some  type  of  decision  logic.6  By  categorizing  the  tracks  into  different  Association  Groups  (AG’s), 
the  number  of  candidate  tracks  to  be  considered  for  association  with  new  measurements  is  reduced,  thus  reducing 
the  computational  load  similar  to  track  partitioning.  This  is  in  essence  a  form  of  coarse  gating  that,  while  useful,  is 
not  all  that  note  worthy.  However,  it  does  lay  the  stage  for  the  second  categorization,  Dwell  Groups  (DG’s).  By  first 
categorizing  the  tracks  into  AG’s,  the  task  of  assigning  membership  to  the  DG  is  simplified.  While  the  computational 
savings  offered  by  AG’s  are  attractive,  the  computational  savings  pale  in  comparison  to  the  potential  savings  of  radar 
resources  offered  by  effective  use  of  DG’s. 

The  savings  with  DG’s  is  seen  when  targets  are  sufficiently  close  that  a  single  track  dwell  is  capable  of  illuminating 
them  so  that  fewer  track  dwells  are  required  to  maintain  track.  This  in  combination  with  measurements  of  opportunity 
from  the  search  operations,  track  on  search,  offers  significant  radar  resource  savings.  A  full  implementation  of  these 
strategies  requires  a  sophisticated  track  management  scheme  that  is  capable  of  organizing  the  tracks  into  the  different 
categories  and  making  efficient  use  of  the  available  measurements.  A  key  feature  of  such  a  tracker  is  the  ability 
to  correctly  suppress  track  dwell  requests  that  have  been  earlier  placed  on  the  scheduling  queue.  These  dwells  are 
considered  redundant  before  the  dwell  requests  are  granted  due  to  measurements  from  other  dwells.  Many  algorithmic 
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solutions  of  varying  sophistication  that  exist  solve  the  fundamental  track  goals.  In  terms  of  radar  resource  efficiency, 
the  better  ones  are  nontrivial. 

The  proposed  approach  offers  notable  improvements  over  simpler  track  management  techniques  when  targets  are 
closcly-spaced.  When  multiple  targets  generate  multiple  detections  from  a  single  dwell,  correctly  associating  them 
with  current  track  files  and  making  full  use  of  the  data  at  hand  is  a  non-trivial  task.  Earlier  approaches  simplified 
this  association  problem  by  considering  only  the  track  that  initiated  the  dwell.  The  closest  measurement  to  this 
track  is  then  selected  and  all  other  measurements  discarded.  Not  only  is  this  an  inefficient  method  in  that  it  does 
not  make  full  use  of  all  measurements,  but  it  is  prone  to  false  tracks  and  cannot  initiate  new  tracks  except  through 
search. 

The  paper  is  organized  as  follows:  Section  two  describes  the  nomenclature  and  assumptions  used  to  describe 
the  problem.  Issues  such  as  the  underlying  assumptions  behind  the  radar  operation  are  discussed  followed  by  the 
definitions  of  the  Association  Group  and  Dwell  Group  as  well  as  track  on  search  (ToS).  Section  three  focuses  on 
Dwell  Group  management  issues  such  as  the  logic  needed  to  properly  implement  this  technique.  The  fourth  section 
presents  some  computer  simulations  to  illustrate  the  improvement  of  the  proposed  technique.  Conclusions  and  areas 
for  future  work  are  discussed  briefly  in  Section  5. 

2.  NOMENCLATURE/DEFINITIONS 

In  order  to  proceed,  it  is  necessary  to  define  the  terms  and  assumptions  pertaining  to  the  problem  at  hand.  This 
is  done  first  with  a  brief  overview  of  the  assumptions  concerning  the  radar  followed  by  descriptions  of  association 
groups,  dwell  groups,  and  track  on  search. 

2.1.  Radar  Assumptions 

The  radar  model  in  question  is  a  notational  multi-faced  phased  array.  The  multiple  arrays  are  knitted  together  into 
a  single  radar  that  is  used  to  perform  search,  acquisition,  and  track.  For  all  three  modes,  the  radar  obtains  range 
and  angle  (azimuth  and  elevation)  measurements  for  each  detection. 


2.2.  Association  Group 

The  definition  of  Association  Groups  is  based  on  the  distance  between  the  tracks  at  a  given  time.  Because  the 
tracks  are  dynamic,  the  membership  in  the  AG’s  must  be  reevaluated  periodically.  The  optimal  schedule  for  such 
reevaluation  is  currently  an  open  issue.  For  simplicity,  it  can  be  assumed  that  this  reevaluation  occurs  on  each  track 
dwell. 

The  distance  used  to  define  the  AG’s  may  be  described  either  geometrically  or  statistically.  If  geometric  distances 
are  to  be  used,  then  the  membership  criteria  is  based  on  arr  sy-face  coordinates.  In  addition  to  the  natural  relationship 
of  array-face  coordinates  to  radar  measurements,  there  is  the  understanding  that  the  targets  must  be  within  at 
least  some  multiple  of  a  beamwidth  of  each  other.  Because  the  beamwidth  in  sine-space  is  independent  of  steer- 
direction,  array-face  coordinates  are  a  convenient  choice.7  The  array-face  coordinates  of  the  track  are  given  as 
tkAFC  =  [r, kx,ky]T  where  r  denotes  range  and  kx  and  ky  denote  the  sine-space  azimuth  and  elevation  angles, 
respectively. 

Let  the  distance  between  any  two  tracks  be  normalized  to  be  either  a  logical  1  if  near  or  a  logical  0  if  far.  Thus, 
given  M  tracks  the  matrix  of  logical  distance,  D  €  fRA/XjV/,  is  given  as 


[D  ];,*  = 


1,  if  | rk  -  Tj |  <  rA  and  -  kVj)2  < 

0,  if  | rk  -  Vj |  >  rA  or  -  kTj)-  +  {kyii  -  kyi  )2  >  0A 


where  ta  denotes  a  range  threshold,  9\  denotes  angle  separation  threshold  in  sine-space,  and  [•]*, j  the  kth,  jth 
component  of  the  matrix,  D.  A  good  choice  for  (9 a  is  the  6  dB  beamwidth  and  the  track  dwell  range  extent  for  r&.  If 
statistical  distances  are  to  be  used,  then  the  logical  distance  is  defined  in  terms  of  track  coordinates  and  covariance 
matrices  in  the  East-North-Up  (ENU)  frame  as 


if  (tkENUk  —  tkENUj)T(PA:  +  P  J  )  1  (tkENUk  —  tkENUj)  <  A 


[Dk*  = 


(2) 


Figure  1.  Association  Group  example. 


where  A  is  a  distance  threshold. 

Once  the  matrix  D  is  found,  the  AG’s  may  be  formed  into  either  non-overlapping  partitions,  or  overlapping  groups. 
The  advantage  of  forming  non-overlapping  partitions  is  easier  implementation.  While  overlapping  groups,  form  more 
AG’s,  they  have  fewer  members  and  thus  require  a  smaller  processing  load  to  evaluate  the  measurement/track 
assignment. 

Intuitively,  the  groups  describe  a  clustering  of  tracks  such  that  there  exists  a  connected  covering  of  AG  members 
with  epsilon  balls  of  radius  defined  by  the  logical  distance  metric.  Thus,  each  AG  contains  members  that  are  potential 
candidates  for  assignment  with  measurements  derived  from  a  dwell  generated  by  any  member  of  the  AG. 

Alternatively,  the  tracks  may  be  partitioned  into  non-overlapping  AG’s  by  considering  the  structure  of  D.  With 
the  proper  bookkeeping,  the  rows  and  columns  of  D  may  be  shuffled  such  that  a  modified  matrix,  D  has  block 
diagonal  form  with  the  blocks  corresponding  to  the  different  AG’s.  A  track  has  membership  into  an  AG  if  it  also  is 
an  index  into  the  corresponding  block  of  D  (see  the  example  described  in  Figure  1).  Membership  into  overlapping 
groups  requires  that  all  tracks  are  within  the  threshold  distance  of  each  other.  Tracks  tkk  and  tkj  are  contained  in 
the  same  overlapping  AG  if  [D ]kj  —  1. 

An  example  is  given  in  Figure  1  where  the  relative  track  positions  are  illustrated  notationally.  Using  either  of  the 
two  distance  measures  discussed  above,  it  is  possible  to  define  the  logical  distance  matrix  shown  below.  Also  shown 
is  the  block  matrix  D.  Both  matrices  are  bordered  by  rows  and  column  of  numbers  that  indicate  the  related  track 
index. 
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Thus,  if  the  AG  are  allowed  to  overlap,  the  following  groups  may  be  formed:  AG\  =  {tk\,  tk2,  tk^:  i/ce } ,  AG 2  =  {^3}, 
and  AG 3  =  {^4, ££5,  tk$}.  From  the  block  structure  of  D,  it  can  be  easily  seen  that  there  are  only  two  groups  that 
do  not  overlap:  AG\  —  {tk\,tk2,tk4,tk$,tk§}  and  AG 2  —  {tk^}. 


2.3.  Dwell  Group 

Dwell  Groups  are  an  attempt  to  make  better  use  of  track  dwells  such  that  fewer  are  needed.  The  concept  is  simple- 
use  a  single  dwell  to  illuminate  multiple  closely-spaced  targets.  Given  a  set  of  track  states  updated  to  the  same  time, 
the  DG’s  are  assigned  such  that  there  exists  a  covering  of  all  tracks  with  dwell  beams.  The  coverage  can  be  defined 
in  terms  of  the  3  dB  beam  contour,  or  a  multiple. 

There  are  two  key  features  of  a  DG-  its  membership  and  its  pointing  direction.  The  membership  is  a  critical 
part  of  the  track  management  logic  that  deals  with  issues  such  as  miss  detections  and  drop/coast  logic.  Once  the 
membership  has  been  decided,  the  actual  beam  pointing  direction  must  be  determined.  While  the  optimal  solution 
is  beyond  the  scope  of  this  paper,  it  must  take  into  account  the  probability  of  detections,  estimated  RCS,  target 
model  dynamics,  and  track  preferences.  One  suboptimal  method  that  is  easy  to  implement  is  to  simply  point  the 
beam  in  the  direction  given  by  the  centroid  of  the  DG  track  members. 
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Figure  2.  Example  of  two  Dwell  Group  memberships  for  the  same  cluster  of  tracks.  The  open  large  circle  represents 
a  dwell  beam  where  the  small  closed  circles  denote  track  positions  and  the  cross-hair  denotes  beam  center. 


Figure  2  illustrates  some  of  the  concerns  in  optimally  assigning  membership  into  the  DG’s  and  in  pointing  the 
dwell  beam  as  well.  In  case  1,  the  memberships  are  not  a  partition  since  tk4  and  tk6  are  in  both  DG1  and  DG2. 
Case  2  is  a  strict  partition.  However,  it  could  be  argued  that  it  is  not  desirable  since  tk4  and  tk6  are  close  to  the 
beam  shoulder.  Case  1,  on  the  other  hand,  does  not  put  these  tracks  much  closer  to  the  beam  center,  but  it  may 
be  possible  to  combine  the  measurements  from  the  two  dwells,  DG1  and  DG2,  to  improve  the  quality  of  the  overall 
measurements. 

Even  though  the  positioning  of  the  DG  may  require  that  some  tracks  fall  into  more  than  one  DG,  it  is  possible 
to  force  a  partition  structure  to  simplify  the  dwell  logic.  This  is  simply  done  by  arbitrarily  assigning  membership  to 
one  DG. 

Given  that  the  optimal  solution  is  currently  beyond  reach,  consider  a  simple  strategy  to  form  the  dwell  groups  and 
position  the  resulting  dwell  beams.  Using  the  same  track  positions  in  Figure  2,  the  simple  algorithm  is  illustrated  in 
three  steps  in  Figure  3.  To  simplify  the  process,  only  tracks  in  an  AG  are  candidates  for  membership  into  the  DG’s 
in  question.  In  other  words,  since  tracks  in  separate  AG’s  are  unlikely  to  associate  with  each  other,  it  is  unlikely  that 
they  are  close  enough  to  be  adequately  illuminated  with  a  few  closely-spaced  beams.  Note  that  the  algorithm  below 
speaks  of  pointing  beams,  it  is  understood  that  no  actual  beams  are  scheduled  while  assigning  membership  to  DG’s. 

step  1:  Form  cluster  of  tracks  of  interest 

Point  beam  at  centroid  of  cluster-  the  tracks  that  are  covered  are  in  DGk 
step  2:  Point  DGk+i  at  the  track  that  is  furtherest  from  DGk 
The  tracks  that  DGk+i  covers  are  in  DGk+i 
Redirect  DGk  beam  at  the  tracks  not  in  DGk+i 
step  3:  Center  the  beams  about  the  DG  members. 

If  all  tracks  are  covered,  stop,  else  eliminate  the  tracks  in  DGk+i  and  go  to  step  1. 

2.4.  Track  on  Search 

Track  on  search  makes  use  of  measurements  collected  during  the  normal  search  operations  of  the  radar  to  update 
existing  tracks.  This  is  made  possible  because  both  range  and  AoA  estimates  are  available  in  search  mode.  While 
it  is  a  good  idea  to  make  use  of  all  available  measurements  from  the  radar,  ToS  measurements  are  not  without 
drawbacks.  In  addition  to  the  added  complexity  in  the  track  management  logic  to  handle  ToS  measurements,  they 
may  also  be  of  reduced  quality  compared  to  normal  track  dwell  measurements.  This  reduced  quality  is  because  the 
search  beam  pointing  location  is  not  directed  at  a  DG.  It  is  not  expected  that  all  members  of  the  DG  will  be  centered 
in  the  beam  such  that  they  all  will  be  sufficiently  illuminated  to  produce  detections.  Similarly,  since  the  quality  of 
many  angle  of  arrival  (AoA)  estimates  degrade  as  the  target  moves  off  beam  center  into  the  shoulder,  the  resulting 
ToS  measurements  will  not  be  of  the  same  quality  as  the  normal  track  dwells.  This  reduction  in  ToS  measurement 
quality  over  track  dwell  measurements  has  not  been  adequately  analyzed  for  the  problem  at  hand  and  may  prove  to 
be  insignificant.  It  is  mentioned  for  completeness. 
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Figure  3.  A  simple  method  of  generating  DG  and  beam  pointing  directions. 

3.  DWELL  GROUP  MANAGEMENT  ISSUES 

As  should  be  expected,  the  improved  efficiency  of  dwell  groups  comes  at  the  expense  of  additional  computational 
complexity.  There  are  more  bookkeeping  issues  in  dealing  with  missed  detections,  coasting,  and  drop  decisions. 
Additionally,  the  adaptive  revisit  time  logic  must  be  reassessed  because  of  the  dynamic  nature  of  DG’s. 

3.1.  Track  Management  Logic 

The  logic  associated  with  this  type  of  track  management  is  a  bit  more  involved  than  with  conventional  methods  in 
which  a  dwell  is  keyed  to  a  track  index.  In  traditional  track  logic,  if  a  dwell  generates  no  measurements  that  associate 
with  its  parent  track,  then  a  decision  is  made  to  either  coast  or  schedule  another  dwell.  However,  with  dwell  groups, 
it  is  not  that  simple  since  several  tracks  are  keyed  to  a  dwell  and  the  membership  can  change  between  dwells.  Some 
of  the  important  issues  that  must  be  addressed  before  the  DG  strategy  can  be  correctly  implemented  are: 

•  acquisition:  the  initiation  of  multiple  new  tracks  within  a  DG  may  require  a  different  dwell  sequence  than 
singletons. 

•  missed  detections:  the  probability  of  miss  detections  of  multiple  tracks  in  a  DG  must  be  understood  as  it 
impacts  coasting  and  drop  track  logic. 

•  coasting  and  drop  track  logic:  decisions  to  drop  track  may  be  influenced  if  the  track  is  a  member  of 
multiple  DG’s.  It  may  be  more  appropriate  to  coast  longer  than  if  the  track  were  a  singleton.  If  the  track  were 
a  singleton,  additional  radar  resources  may  be  needed  to  re- acquire  than  if  it  were  coasted  in  a  DG. 

•  merge  track  decisions:  when  to  join  multiple  tracks  together  should  depend  on  DG  membership. 

•  scheduling  track  dwells:  DG’s  influence  the  track  beam  positions  and  their  timing.  To  fully  capture  the 
efficiencies  offered  by  DG’s,  it  is  necessary  to  have  proper  logic  that  will  inhibit  previously  scheduled  dwells 
that  are  considered  redundant  because  of  measurements  from  other  nearby  dwell  groups  that  associate  with 
the  tracks. 

There  exist  many  possible  choices  regarding  the  implementation  details  for  these  items  that  will  result  in  a  working 
tracker,  and  a  detailed  discussion  will  not  be  given  here. 

3.2.  Next  Dwell-time  Request  (Adaptive  Revisit) 

In  scheduling  dwell  times,  there  are  two  choices:  adaptive  revisit  time  and  periodically  scheduled.  In  attempting 
to  make  the  most  economical  use  of  dwells,  adaptive  revisit  time  is  the  better  choice.  The  adaptive  revisit  time  is 
computed  to  maximize  the  time  between  dwells  such  that  the  target  dynamics  model  is  constrained  to  be  illuminated 
by  the  dwell  when  it  is  sent.5  Issues  that  must  be  taken  into  account  include  the  beamwidth,  range  extent  of  the 
dwell,  and  accuracies  of  the  target  state  estimates  as  well  as  the  maximum  thrust  and  turn  acceleration  of  the  target. 

Since  the  targets  are  free  to  maneuver,  the  DG  membership  may  change  before  the  dwell  if  the  adaptive  revisit 
period  ( i.e .  time  between  two  consecutive  DG  dwells)  is  too  long.  Because  of  poor  track  maintenance  performance 


observed  through  simulation  studies,  it  was  noted  that  this  should  be  avoided.  Thus,  if  it  is  observed  that  the  DG 
memberships  change  before  the  next  scheduled  dwell,  the  adaptive  time  is  reduced  by  roughly  half. 

To  further  simplify  the  tracker’s  task,  adaptive  revisits  are  only  allowed  for  DG’s  with  only  one  member  ( i.e . 
singleton).  These  singletons  can  capitalize  on  the  throughput  gains  offered  by  adaptive  revisit  times.  For  DG’s  with 
more  than  one  member,  periodically  scheduled  revisit  times  proved  to  provide  more  robust  track  maintenance. 

4.  COMPUTER  SIMULATIONS 

To  illustrate  the  concept  of  the  proposed  AG/DG  strategy,  computer  simulations  are  presented.  The  radar  model 
used  to  derive  the  plots  is  a  notational  ship  based  radar  of  medium  fidelity  level  .  The  radar  operates  at  3.4  to  4.0 
GHz  with  four  planar  arrays  that  provide  360°  of  azimuthal  coverage  and  75°  in  elevation.  Each  of  the  four  arrays  is 
tilted  back  15°  and  faces  90°  apart  from  each  other  to  give  full  azimuthal  coverage.  The  field  of  view  of  each  array  is 
such  that  any  two  arrays  overlap  in  azimuth  by  about  10°.  Linear  motion  of  the  platform  is  modeled  as  well  as  roll, 
pitch,  and  yaw.  The  four  arrays  are  knitted  together  into  a  single  radar  that  is  used  to  perform  search,  acquisition, 
and  track.  For  simplicity  the  same  suite  of  waveforms  are  assumed  used  for  each  task  with  individual  waveforms 
varying  in  frequency  and  pulsewidth.  After  a  track  has  been  established,  the  proper  pulsewidth  is  selected  to  adjust 
the  energy  on  target.  The  range  estimates  are  based  on  a  discrete-phase  coded  pulse  compression  model  while  the 
AoA  measurements  are  generated  with  sum  and  difference  monopulse  for  both  azimuth  and  elevation.  As  a  result, 
if  range  resolution  is  lost,  the  monopulse  estimates  will  be  corrupted  due  to  the  fact  that  multiple  targets  within  a 
single  range  cell  contaminate  monopulse  AoA  estimates.  Monopulse  estimates  are  available  in  all  modes  of  operation 
such  that  data  collected  during  search  may  be  efficiently  used  in  ToS. 

High  fidelity  issues  such  as  element  level  effects  are  not  modeled,  instead  an  analytic  expression  is  used  to  represent 
the  sum  and  difference  antenna  patterns.  This  level  of  fidelity  allows  issues  such  as  multiple  detections  in  adjacent 
search  beams  to  be  studied  along  with  attenuation  effects  due  to  off-boresight  beams. 

Figures  4-6  present  different  views  of  the  same  simulation  run.  The  simulation  starts  with  two  aircraft  in  flight, 
with  the  ship  denoted  by  the  S  at  location  (0,0)  as  seen  in  Figure  4.  The  aircraft  then  join  and  fly  in  formation  until 
the  end  of  the  simulation.  In  ail  of  these  figures,  the  truth  trajectories  are  denoted  with  a  broken  line  while  the  solid 
lines  represent  tracks. 

Figures  5  and  6  depict  sub-regions  of  Figure  4  to  illustrate  finer  detail.  These  figures  also  show  beam  pointing 
track  dwells  and  the  resulting  measurements  along  with  ToS  measurements.  The  ToS  dwells  are  not  shown,  because 
there  is  no  equivalent  range  component  of  the  search  dwell-  it  does  not  carry  the  same  meaning  as  it  does  in  the 
track  dwell  and  would  clutter  the  figures.  The  track  dwell  beam  pointing  directions  for  tracks  1  and  2  are  shown 
with  circles  and  diamonds  while  x  and  +  denote  measurements  and  the  square  is  used  to  show  ToS  measurements. 

In  Figure  5  the  targets  are  in  separate  AG’s  and  separate  DG’s  resulting  in  separate  track  dwells  for  each  target. 
Adaptive  revisit  time  is  used  so  that  the  spacing  of  the  track  dwells  is  aperiodic.  The  ToS  dwells  which  are  not 
shown  are  periodic.  Their  resulting  measurements  may  be  seen  to  overlap  the  track  dwells  on  occasion.  It  is  these 
times  where  the  ToS  dwells  may  be  considered  redundant  to  the  track  dwells,  because  the  search  dwells  must  be 
scheduled  and  never  inhibited.  Conversely,  isolated  ToS  measurements  may  be  seen  where  the  search  dwell  happened 
just  before  a  scheduled  track  dwell,  inhibiting  it  and  giving  some  of  the  promised  radar  resource  efficiency.  Also  of 
note  is  where  the  targets  transition  from  one  array  face  to  another.  This  array  face  hand-off  can  be  seen  as  the  two 
closely-spaced  ToS  measurements  on  each  track. 

A  closeup  of  the  region  where  the  aircraft  come  together  is  shown  in  Figure  6.  The  targets  are  close  enough  to 
be  in  the  same  AG  since  Track  1  dwells  are  suppressed  because  Track  2  is  being  fed  by  measurements  from  Track  1 
dwells.  The  two  tracks  are  in  the  same  DG,  DGi,  at  the  region  where  the  Track  1  dwells  straddle  the  two  tracks. 
Because  of  an  artifact  due  to  the  plot’s  scale,  it  appears  that  the  dwell  favors  Track  2,  but  it  really  is  in  the  middle. 
Also  note  that  the  track  logic  implemented  in  the  simulation  moves  from  adaptive  to  periodic  revisits  when  the  DG 
has  more  than  one  member. 


5.  CONCLUSIONS /FUTURE  WORK 

A  new  track  management  technique  that  makes  efficient  use  of  radar  resources  has  been  presented  along  with 
computer  simulations  to  illustrate  its  operation.  Key  features  of  this  strategy  involve  the  use  of  association  groups 
and  dwell  groups.  Several  different  methods  of  structuring  the  AG’s  and  DG’s  were  presented  along  with  track  logic 
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Figure  4.  Two  tracks  merging  with  AG/DG  tracker  and  ToS  enabled,  AG  uses  geometric  distance. 


Track  =  AG/DG,  Revisit  =  adaptive  Assign  =  MJV-pure  ToS  =  on  AG  =geometric 


13 


14  15  16 

East  (km) 


19 


20 


Figure  5.  Two  tracks  merging  with  AG/DG  tracker  and  ToS  enabled,  AG  uses  geometric  distance. 
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Figure  6.  Close  up  view  of  tracks  after  the  targets  join  and  the  AG/DG  technique  comes  into  play. 

and  management  issues  that  are  of  concern  when  implementing  these  methods.  Many  issues  are  left  open  and  promise 
to  be  rich  areas  of  future  research.  Performance  analysis  of  the  different  AG/DG  structures  discussed  along  with  the 
optimal  DG  beam  positions  are  obvious  areas  for  further  study.  Also  a  detailed  performance  analysis  is  needed  to 
compare  the  new  algorithm  with  existing  ones  with  regard  to  track  metrics  and,  of  course,  radar  resource  efficiency. 
Other  areas  of  interest  involve  further  refinement  of  the  management  logic  to  make  the  tracker  with  respect  to  radar 
resources  more  efficient  and  robust  to  both  noise  and  smart  jammers. 
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Extended  Abstract 

Monopulse  is  a  simultaneous  lobing  technique  for  determining  the  angular  location  of  a 
source  of  radiation  or  of  a  “target”  that  reflects  part  of  the  energy  incident  upon  it  [1],  In  an 
amplitude  comparison  monopulse  radar  system,  a  pulse  of  energy  is  transmitted  directly 
at  the  predicted  position  of  the  target,  and  the  target  echo  is  received  with  two  beams 
that  are  squinted  relative  to  the  predicted  position  of  the  target  (note  that  two  beams  are 
required  for  each  angular  coordinate).  Traditionally,  the  Direction-Of- Arrival  (DOA)  of  a 
target  is  estimated  with  the  in-phase  part  (i.e.,  the  real  part)  of  the  monopulse  ratio,  which 
is  formed  by  dividing  the  difference  of  the  two  received  signals  by  their  sum.  However,  the 
presence  of  reflections  from  the  sea  surface  in  the  main  lobe  of  the  antenna  pattern  can 
cause  severe  errors,  when  using  the  monopulse  ratio  as  the  target  DOA  estimate.  While  the 
problem  of  tracking  low  elevation  targets  has  been  studied  extensively  in  the  last  30  years 
[1-8],  the  statistics  of  monopulse  measurements  of  a  target  in  the  presence  of  sea-surface 
induced  multipath  were  not  available  in  the  literature  until  [9,10]. 

In  tracking  of  low  elevation  targets  with  a  monopulse  radar,  the  presence  of  reflections 
from  the  sea  surface  can  cause  severe  errors  in  the  DOA  measurements  of  the  target. 
Since  the  target  echoes  that  are  received  directly  and  via  the  sea  surface  are  unresolved 
in  time  and  frequency,  tracking  targets  in  the  presence  of  sea-surface  induced  multipath 
is  a  special  case  of  tracking  unresolved  targets.  The  sea-surface  reflection  is  modeled  by 
a  specular  (coherent)  component  and  a  diffuse  (noncoherent)  component.  The  probability 
density  function  (PDF)  of  the  measured  amplitude  of  the  sum  signal  and  the  amplitude- 
conditioned  means  and  variances  of  the  in-phase  and  quadrature  monopulse  ratios  are  given 
in  [9,10]  for  low  elevation,  fixed-amplitude  targets  in  the  presence  of  sea-surface  induced 
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multipath-.  The  PDF  and  the  means  and  variances  include  the  effects  of  both  diffuse  and 
specular  multipath  and  the  means  and  variances  of  the  monopulse  ratios  are  studied  for 
radars  operating  at  S,  X,  and  Ku  bands. 

The  diffuse  reflection  coefficients  in  [9.10]  that  were  taken  from  [11]  were  found  to  be  in¬ 
consistent  with  the  diffuse  reflection  coefficients  given  by  [12].  The  use  of  diffuse  reflection 
coefficients  of  [11]  has  been  found  to  exaggrate  the  impacts  of  the  diffuse  multipath  com¬ 
pared  to  the  effects  found  in  real  radar  data.  In  this  paper,  the  inconsistencies  in  the  diffuse 
reflection  coefficients  will  be  reconciled  and  the  means  and  variances  of  the  monopulse  ra¬ 
tios  will  be  studies  for  various  operating  frequencies.  The  discrepancies  between  predicted 
and  actual  monopulse  measurements  for  targets  in  the  presence  of  multipath  in  [4]  are 
thought  to  be  the  result  of  diffuse  multipath,  and  the  paper  include  the  results  of  an 
investigation  into  these  discrepancies. 

Additional  Background 

The  sea-surface  reflection  is  modeled  with  a  specular  (coherent)  component  and  a  diffuse 
(noncoherent)  component.  The  specular  reflection  is  caused  by  a  smooth  (‘‘mirror- like") 
surface  and  the  diffuse  reflection  is  caused  by  the  surface  irregularities.  While  the  spec¬ 
ular  reflection  coefficient  is  a  deterministic  number  which  depends  on  several  unknown 
parameters,  the  diffuse  reflection  has  a  random  nature  that  is  often  modeled  as  a  complex 
Gaussian  process.  Generally,  the  sea  surface  is  perturbed  by  small  irregularities,  and  both 
reflection  components  are  present. 

The  signal  received  from  a  low  elevation  target  in  the  presence  of  sea-surface  induced 
multipath  includes  four  components.  The  first  part  travels  directly  to  the  target  and 
returns  directly  to  the  radar,  while  the  second  part  travels  to  the  target  via  the  sea  surface 
and  returns  directly  to  the  radar.  The  third  part  travels  directly  to  the  target  and  returns 
to  the  radar  via  the  sea  surface,  while  the  fourth  part  travels  to  the  target  via  the  sea 
surface  and  returns  to  the  radar  via  the  sea  surface.  In  the  presence  of  sea-surface  induced 
multipath,  the  in-phase  and  quadrature  portions  of  the  sum  and  difference  signals  for 
elevation  measurements  are  given  by 

sj  =  at  cos (j>  +  2 atgps  cos {<f>  +  A<j>)  +  at(gps)2  cos (cf>  +  2A 4>)  +  2atgpd  cos (<fid  +  A<j>) 

+  at{gpd)2  cos(24>d  +  2A(f)) +  nSI  (1) 

sq  =  at  shuf)  +  2 atgps  sin (<t>  +  A <f>)  +  at{gps )2  sin (<j>  +  2A</>)  +  2 atgpd  sin(<^  +  A</>) 
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(2) 


+  ott{gpd )2  sin(2 (f>d  +  2 A <f>)  +  nSQ 

di  =  atih  cos <j>  +  at(rjt  +  rji)gps  cos (<f>  +  A<f>)  +  atpi{gps)2  cos (<j>  +  2A</>) 

+  Mvt  +  Vi)gpd  COS {<t>d  +  A(p)  +  atru (gpd)2  cos(20d  +  2A<£)  +  ndI  (3) 

dQ  =  OitT}t  sin</>  +  Oit{rjt  +  Vi)gps  sin (<f>  +  A <j>)  +  atpi{gps)2  sin (<j>  +  2A0) 

+  at(rjt  +  T)j)gpd  sin (</>d  +  A<f>)  +  atgi{gpd)2  sin(2^d  +  2A <f>)  +  ndQ  (4) 

where 


at  —  \/KAtG^(6t)po 

k  =  proportional  to  the  transmitted  power 
At  =  voltage  amplitude  of  the  target  backscatter 
Gs(9)  =  sum  channel  antenna  gain  at  angle  6 
Ga  (9)  =  difference  channel  antenna  gain  at  angle  9 
9t  =  off-boresight  angle  of  the  target 
Po  =  matched  filter  gain 


4>  =  phase  of  the  directly  returned  signal  echo 
r]t  =  tt4^4  =  DOA  of  the  target 

Gz( Pt) 

Pi  ~  77^777^  ~  DO  A  of  the  target’s  image 
G  y,[9i) 

9 i  =  off-boresight  angle  of  the  target’s  image 
_  Gz(9/) 

9~  Gx(9t) 

A<f>  =  phase  difference  between  the  direct  and  specular  reflections 


DOA  of  the  target’s  image 


4>d  =  uniformly  distributed  phase  of  the  diffuse  reflection 


ps  =  specular  reflection  coefficient 


Pd  =  Rayleigh  diffuse  reflection  coefficient  with  parameter  pd o 
nsi  ~  N{ 0,  erf)  nSQ  ~  N(0,  erf) 

ndi  ~  N{ 0,  erf)  ndQ  ~  N( 0,  <rd) 


The  at  is  the  voltage  amplitude  of  the  target  echo  in  the  absence  of  multipath  reflections. 
Note  that  A <j>  includes  the  phase  difference  due  to  both  the  Path  Length  Difference  (PLD) 
and  the  specular  reflection  at  the  sea  surface,  which  is  approximately  n.  The  receiver  errors 
nsi,  nSQindii  and  ndQ  are  assumed  independent.  Both  ps  and  pdo  depend  on  the  sea  state, 
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properties  of  the  seawater,  polarization  of  the  transmitted  waveform,  grazing  angle  at  the 
point  of  the  sea-surface  reflection,  and  wavelength  A  of  the  carrier. 

The  first  term  on  the  right  side  of  (1)  and  (2)  corresponds  to  the  echo  received  directly 
from  the  target,  while  the  second  and  third  terms  correspond  to  the  three  echoes  that 
are  the  result  of  the  specular  reflection  at  the  sea  surface.  The  fourth  and  fifth  terms 
represent  three  echoes  that  result  from  the  diffuse  reflections  at  the  sea  surface.  The  phase 
information  associated  with  (p  in  the  fourth  and  fifth  terms  is  lost  due  to  the  presence  of 
the  random  phase,  (pd-  In  the  difference  signals  of  (3)  and  (4),  the  second  and  fourth  terms 
on  the  right  side  of  (3)  and  (4)  include  echoes  from  two  different  DOAs,  r)t  and  tjj. 


The  specular  reflection  coefficient  is  computed  as 


where 


9o  — 


r  = 


rH 

Tv 


(  T  exp(-87 r2^),  0  <  g0  <  0.1 

Ps  = 

{  „  0.81254 

1  1  +  8^5' 

(5) 

CTh  . 
—  sin 
A 

'tpga 

(6) 

( rv, 

vertical  polarization 

\th, 

horizontal  polarization 

sin  ipga  -  ■v/ec  -  cos 2i\,a 

sin  Ipga  +  \ftc~  COS2 'Ip 


ga 


ec  sin  ipga  -  sjjc  ~  cos2^. 


9“ 


— 


ec  sin  xpga  +  y/tc~  cos 2t/>aa 
es-to  .{(es~eo)uc T  _  4?T£M 


1  +  U)2T2 


+  u£t 


2t-2 


uc  —  —  =  carrier  frequency 
A 

<?h 
Ipga 

es 


wavelength  of  the  carrier  frequency 

RMS  sea-surface  elevation  above  the  mean  level 

grazing  angle 

static  dielectric  parameter  of  the  seawater 
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eo  =  4.9  for  seawater 
t  =  relaxation  time  of  the  seawater 
(7j  =  ionic  conductivity  of  the  seawater 


According  to  [11],  the  Rayleigh  parameter  for  the  diffuse  reflection  coefficient  is  computed 


\/2|r|3.68<fc,, 


0  <  gQ  <  0.1 


Pdo  =  {  y/2\r\ (0.454  -  0.858&,),  0.1  <  ga  <  0.5 


v/2|r|0.025, 


9o  >  0.5 


Letting  A  and  ip  denote  the  measured  amplitude  and  phase  of  the  sum  signal  gives 


sj  =  A  cos  ip  sq  —  A  sin  ip 


Also,  let 


2  <j| 


where  K  denotes  the  Signal-to-Noise  Ratio  (SNR)  of  a  fixed-amplitude  target  in  the  absence 
of  multipath  and  3R0  denotes  the  observed  SNR. 

Let  $  denote  the  parameters  {at,  A<p,  ps,  Pdo,&s,Vd}-  With  E[-]  denoting  the  expected 
value  and  VAR[-]  denoting  the  variance,  then  [9,10] 


£[5ft0|$]  =  [(1  -I-  2 psg  cosAp  +  p2sg 2)2  +  8p2d0g2(l  +  p2d0g2)]U  +  1 
VAR[^|$]  =  [8p2d0g2(l  +  pd0g2)$t  +  1 

x  [ [2(1  +  2 psg  cosA <p  +  p2sg2)2  +  8p2dQg2{\  +  p2d0g2)]$  +  1 


Denoting  s  =  si  +  jsQ  and  d  =  dj  +  jdQ,  the  in-phase  and  quadrature  parts  of  the 
monopulse  ratio  are  given  by 


(ds  diSi  +  SQdQ 

yi  =  M~)=  .2  I  ,2  ~ 

s  SI  '  SQ 

Md,.  dnsi  -  disn 

V  =  j_  „2 


Si  +  Sq 
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The  expected  value  of  the  monopulse  measurements  [9.10]  are  given  by 


E[yi |A,  $]  =  ^  +  /1|0(^(1  +  2 psg  cosA</>  +  p|t/2))  ^  +  (pstf)2^,  -  ~) 

(14) 

(15) 


+  (vt  +  m  ~  2  — )ps5  cosA0 
V  Pn' 

£[i/q|A,$]  =  /i|o(  — (1  +  2 psg  cos  Atfi  +  p2sg2)) 

\  Vn  ' 


atPsg 


A 


(jli  ~  m)  sinA cj) 


where 


P12  2  a2p^0g2[r/t  +  (1  +  2p^2)r;j] 

Pn  4a2p20<72[l  +  p20<?2]  +  t72 


(16) 

(17) 


The  Iq{x )  and  I\(x )  are  the  zero  and  first  order  modified  Bessel  function  of  first  kind. 
The  form  of  (14)  suggests  that  separating  the  bias  in  the  monopulse  ratio  (i.e.,  as  a 
DOA  estimate  of  the  target)  into  diffuse  and  specular  components  as  in  [8]  may  not  be 
appropriate.  For  example,  when  measured  amplitude  of  the  sum  signal  is  rather  small, 
the  bias  in  the  monopulse  ratio  will  result  from  the  diffuse  multipath.  On  the  other  hand, 
when  the  measured  amplitude  of  the  sum  signal  is  very  large,  the  effects  of  the  diffuse 
multipath  will  be  very  small. 


Since  the  statistics  of  the  monopulse  ratios  as  presented  above  are  a  function  of  the  mea¬ 
sured  amplitude,  A,  of  each  pulse  or  subpulse,  the  statistics  will  be  illustrated  by  using 
the  expected  value  of  the  observed  SNR  in  (10)  for  the  5R0.  Using  the  $ft0  =  E[JR0|$]  of 
(10)  and  the  sea-surface  reflection  model  above,  trajectories  of  E[yi |A,  $],  E[yq | A,  $],  and 
the  variances  were  generated  for  low-elevation  targets,  with  the  antenna  boresight  pointing 
directly  at  the  target.  Trajectories  were  generated  for  a  radar  operating  at  4  GHz  in  S 
band.  For  all  cases,  the  sea  surface  was  assumed  to  have  a  RMS  wave  height  of  0.25  m, 
and  the  radar  is  vertically  polarized  and  20  m  above  the  sea  surface.  Also,  for  all  cases, 
the  use  of  Sensitivity  Time  Control  (STC)  will  be  assumed  so  that  the  SNR  of  the  target 
only,  Ift,  is  independent  of  range. 

For  the  S  band  case,  the  radar  was  modeled  to  have  a  one-way  bcamwidth  of  2.5°  and  a 
squint  angle  of  1.05°.  The  target  travels  from  a  range  of  20  km  to  a  range  of  5  km  at  an 
altitude  80  rn.  The  means  of  yi  and  yq  and  the  associated  standard  deviations  are  shown 
in  Figures  1  through  3  for  a  16-dB  target  (i.e.,  16  dB  in  the  absence  of  multipath).  The 
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Figure  1.  Means  of  Monopulse  Ratios  for  16-dB  Target  at  an  Altitude  of  80  m  and  S  Band 


Figure  2  Means  of  the  Complex  Monopulse  Ratios  for  16-dB  Target  at  an  Altitude  of  80 
m  and  S  Band 

oscillatory  form  of  yi  and  \)q  in  Figure  1  is  the  result  of  changes  in  A 4>.  The  dips  in  the 
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0.9 


SOLID:  IN-PHASE  RATIO 
DASH:  QUADRATURE  RATIO 


Figure  3.  Standard  Deviations  of  Monopulse  Ratios  for  16-dB  Target  at  an  Altitude  of  80 
m  and  S  Band 

peaks  of  E[yj |A,<3>]  are  the  result  of  the  diffuse  multipath  dominating  the  return  signal 
when  the  direct  and  reflected  signals  almost  cancel.  The  results  of  Figures  1  are  given  in 
the  complex  plane  in  Figure  2  along  with  the  same  results  for  pdo  —  0.  Note  that  Figure  2 
gives  only  the  results  for  the  target  ranges  of  15  km  to  8  km.  The  dash  lines  in  Figure  2 
give  the  means  of  yi  and  yQ  for  specular-only  reflections  at  the  sea  surface.  Figure  2  shows 
that  ignoring  the  effects  of  the  diffuse  reflections  will  result  in  significant  errors  in  the  DOA 
estimation  for  a  target  in  the  presence  of  sea-surface-induced  multipath.  Note  that  these 
curves  utilize  teh  diffuse  reflection  coefficients  of  [11].  Based  on  the  examination  of  actual 
radar  data,  the  effects  of  the  diffuse  multipath  appears  to  be  too  large.  Figure  3  shows 
that  the  standard  deviations  of  the  in-phase  monopulse  ratios  are  significantly  less  when 
the  target  is  out  of  the  multipath  nulls  (i.e.,  peaks  in  the  expected  value  of  the  in-phase 
monopulse  ratios)  and  at  longer  ranges,  where  the  specular  reflections  are  relatively  strong, 
and  the  angles  between  the  target  and  the  image  are  rather  small. 
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A  bs  tract- 1  he  steady-state  gains  and  error  covariance  are 
derived  for  a  two-state  Kalman  filter  (i.e.,  an  a ,  3  filter) 
for  tracking  with  linear  frequency  modulated  (LFM)  wave¬ 
forms.  A  procedure  is  given  for  calculating  the  tv,  3  gains 
from  the  tracking  index,  T,  the  sample  period,  T,  and  the 
range-Doppler  coupling  coefficient,  AT  The  steady-state 
error  covariance  is  found  to  be  a  simple  function  of  tv,  3, 
A t,  T  and  the  measurement  noise  variance,  cr~.  A  gain 
scheduling  technique  for  tv  and  3  during  filter  initializa¬ 
tion  is  also  given. 


where  Wk  ~  ;V(0 yQk)  is  the  process  noise,  defines  a 
linear  constraint  on  the  dynamics  and  GT  is  the  input 
matrix  for  maneuver.  The  target  state  vector  AT  contains 
the  range  and  range  rate  of  the  target  at  time  k.  The  linear 
measurement  model  is  given  by 

Yk  —  Hk  AT  +  Vk  (2) 

where  Yk  is  typically  the  target  position  measurement  and 
Vk  ~  N(Q,Rk)  is  the  measurement  noise,  and  Hk  is  the 
output  matrix  that  will  include  range-Doppler  coupling 
coefficient. 


I.  INTRODUCTION 

The  Kalman  filter  produces  an  optimal  estimate  of 
the  target  state  when  given  the  motion  model  and  a  se¬ 
quence  of  sensor  measurements  corrupted  with  white  Gaus¬ 
sian  errors.  The  computational  burden  of  maintaining  the 
Kalman  filter  may  prohibit  its  use  when  many  targets  are 
being  tracked.  Steady-state  filters  are  also  good  for  filter 
design  and  analysis.  Constant  data  rate  and  constant  error 
covariances  required  for  steady-state  filters  are  not  typi¬ 
cally  true  in  real  systems,  so  approximate  gains  are  typi¬ 
cally  used.  The  gains  are  based  on  the  steady-state  gains 
of  the  Kalman  filter  and  implemented  either  using  a  fixed 
gain  or  an  easily  generated  gain  schedule.  The  resultant  fil¬ 
ter  for  tracking  range  and  range  rate  is  called  the  tv,  3  filter, 
where  approximate  filter  gains  cv,  3  may  be  used  instead 
of  the  optimal  Kalman  gains  [I].  In  [2],  the  range-Doppler 
coupling  of  linear  frequency  modulated  (LFM)  waveforms 
is  shown  to  have  a  significant  effect  on  tracking  accuracy. 
In  the  Appendix  of  [1]  and  in  [3],  the  steady-state  gains 
and  error  covariance  are  derived  for  range  measurements 
without  range-Doppler  coupling.  In  this  paper,  the  range- 
Doppler  coupling  associated  with  LFM  waveforms  is  in¬ 
cluded  in  the  measurement  equation,  and  expressions  for 
the  steady-state  gains  and  error  covariance  are  calculated. 


The  measured  range  and  the  variance  for  a  LFM 
waveform  and  the  range-Doppler  coupling  coefficient  for 
a  LFM  waveform  are  given  by  [2]  as 


yk  =  rk  +  Ati'k  +  Vk 


=  Etf]  = 

A  /  = 


A  t’c2 

W™ 

for  _  for 

B 


h-fi 


c~ 


(3) 

(4) 


(5) 


where  t‘k  is  the  true  range  at  time  Ar,  R  is  the  signal-to- 
uoise  ratio,  c  is  the  speed  of  light,  /o  is  the  nominal  carrier 
frequency,  f\  is  the  initial  frequency,  and  /•>  is  the  final 
frequency.  Also  \8\  is  the  bandwidth  of  the  LFM  waveform 
and  r  is  the  pulse  length.  Note  that  the  demodulation 
frequency  can  be  modified  with  the  range  rate  estimate  to 
reduce  the  bias  in  the  range  measurement  due  to  range- 
Doppler  coupling.  However,  if  the  range  rate  estimate 
is  used  to  remove  the  bias  in  the  range  measurements, 
the  measurement  errors  become  correlated  with  the  state 
estimate  and  the  typical  Kalman  filter  must  be  modified 
to  account  for  the  state  correlated  measurement  errors. 


The  Kalman  filtering  equations  associated  with  the 
state  model  in  (1)  and  the  measurement  model  in  (2)  are 
given  by  the  following  equations. 


II.  MATHEMATICAL  MODEL 


Time  Update: 


The  dynamics  model  commonly  assumed  for  a  target 
in  track  is  given  by 

AT-+1  =  YkXk  +  Gku.’k  (1) 


'This  research  was  accomplished  through  funding  from  the  Fac¬ 
ulty  Leadership  Program  of  the  Georgia  Tech  Research  Institute, 
Atlanta  Georgia. 


ATife-1  =  Fk-xXk-i\k-x  (h) 

-  Fk-iPk-i\k-i^k~i  +  Gk-iQk-\Gl-\  (7) 

Measurement  Update: 

A*  =  Pk\k-,Hl[HkPk\k_yHl  +  Rk}~1  =  Pk]kH{  R^(S) 

-^k\k  —  +  Afc[y*  —  (9) 

[\\k  =  [I-  l<kHk}Pk\k-i  (10) 


where  the  subscript  notation  {k\j)  denotes  the  state  esti¬ 
mate  for  time  tk  when  given  measurements  through  time 
tj .  The  Kk  is  the  Kalman  gain  that  minimizes  the  mean 
square  error  in  the  state  estimate.  Also,  note  that  Xk  — 
N (Xk\k,Pk\k)  with  Xk\k  and  Pk\k  denoting  the  mean  and 
error  covariance  of  the  state  estimate,  respectively. 

III.  THE  tv,  3  FILTER 


Inserting  (7)  into  (10)  and  setting  Pk\k  =  Pk_i\k_i  =  P 
for  steady-state  conditions  gives 

[/  -  KH]~l  P  =  FPFr  +  GGl'al  (22) 

Equating  the  (1,2)  and  (2.2)  elements  of  (22)  gives 
expressions  for  pyt  and  Substituting  p\2  and  into 
(21)  gives 


For  the  Kalman  filter  in  steady-state  conditions, 
Pk\k  =  Pk-i\k-u  Pk+i\k  =  Pk\k~u  and  Kk  =  Kk - 1 ■  For  a 
Kalman  filter  to  achieve  these  steady-state  conditions,  the 
error  processes,  wk  and  vk ,  must  have  stationary  statistics 
and  the  data  rate  must  be  constant.  When  the  noise  pro¬ 
cesses  are  not  stationary  or  the  data  rate  is  not  constant, 
a  filter  using  the  steady-state  gains  will  provide  subopti- 
rnal  estimates.  The  a,  3  filter  is  the  steady-state  Kalman 
filter  for  tracking  nearly  constant  range  rate  targets.  The 
a,  3  filter  is  a  single  coordinate  filter  that  is  based  on  the 
assumption  that  the  target  is  moving  with  constant  range 
rate  plus  zero-mean,  white  Gaussian  acceleration  errors. 
Given  this  assumption,  the  filter  gains  cv  and  3  are  cho¬ 
sen  as  the  steady-state  Kalman  gains  that  minimize  the 
mean-square  error  in  the  range  and  range  rate  estimates. 

For  the  cv,  3  filter  for  piecewise  constant  acceleration 
errors. 


A'k  -  [r*;  ffc 

1 


1 T 


Fk  -  F  — 

Gk  -  G  - 
Hk  =  H  = 


T 
0  1 

X* 

2 


2’] 


[1  At] 
Rk  =  R  =  a'i 
Qk  =  Q-  al 

"1T 


Kk 


=  K  =  [  a  - 
L  t 


(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 


where  rk  and  rk  are  the  range  and  range  rate  of  the  target, 
respectively,  a2  and  are  the  variances  of  measurement 
noise  and  process  noise,  respectively,  T  is  the  sample  pe¬ 
riod  between  measurements,  and  At  is  the  range-Doppler 
coupling  coefficient. 


Let  the  steady-state  error  covariance  matrix  of  the 
filtered  estimates  for  the  cv,  3  filter  be  denoted  as 


Pk+i\k+i  —  Pk\k  —  P  — 


P  ii 
Pl2 


V 12 
P22 


(18) 


Using  (14),  (15),  and  (18)  in  (8)  gives  the  steady-state  gain 
as 


Kk  =  Kk-i  =  K  =  PHrR~l  = 


Pn  +Pil>A< 
Pv>  +  P'22^t 


1  (19) 


Using  (17)  with  (19)  gives 

«  =  (Pn  +  pi:>A*)cr“2 

3  —  [pn  +  P'>‘>&t)Tcr‘2 


r 


3'2 


1  —  cv  —  3  — 

T 


(26) 


where  F  is  the  tracking  or  maneuvering  index.  Using  (26) 
to  eliminate  cr2,  in  the  expressions  of  py>  and  gives 


PV2 


Pn  = 


3*1 


fal[2a-0) 

2 T*  (l-a-  0^) 


Substituting  py>  into  (20)  gives 

■>  3&i 


V 11  =  “ 


1  L 


2cv  T  3  At 

~~ ~ )  T 


(27) 


(28) 


1  —  CV 


2cv  -|-  3  \  Xt 

~  2  )  T 


A  t 

Y 

(29) 


The  steady-state  error  covariance  is  then  given  by  (27) 
through  (29).  Equating  the  (1,1)  elements  of  (22)  and 
using  (27)  through  (29)  gives 

ti's3^  ~F  u>2$*  T  ai3  -h  no  —  0  ('-19) 

where 


{At\- 

< 3 

4  T 

(31) 

„  f  AG 

v  "  A  t  A  t  1  —  cv 

“ 3  =  “2a  [t, 

)  -2“r+8-r+  4 

(32) 

•>  At 

At 

“l  =  -3a2—  - 

■  a“  +  2a  —  +  3a  —  2 

(33) 

2  3 

a  u  =  a  —  a 

(34) 

There  are  three  solutions  for  (30).  By  comparing  with  the 
Kalman  filter  numerical  results,  the  valid  solutions  of  (31) 
for  3  are  given  by 


(20) 

(21) 


x 


(35) 


For  —  >  0,  both  solutions  of  (35)  are  valid  for  3 

At 

depending  on  the  selection  of  F  and  — .  The  (i.e., 

minus  sign  in  the  front  of  the  radical)  is  used  for  3  with 
increasing  cv  until  the  maximum  cv  is  reached,  which  is 
given  by 


Wmax  — 


2A  t 
~T 


+  1 


4  A  t 


+  1 


2  |  — 
i  T 


(36) 
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Then  3+  (i.e.,  plus  sign  in  the  front  of  the  radical)  is  used 
for  3  from  cvmar  with  decreasing  «  and  it  is  valid  as  long 

as  1  —  q  —  >  0,  that  is,  F  in  (26)  is  a  real  number. 

For  ^  <  0,  ,3-  is  the  only  solution  for  0.  If  ^  —  —  i, 
(30)  becomes  a  quadratic  equation  given  by 


b202  +  bx0  +  60  =  0 

(37) 

cv  1 

62=8-2 

(38) 

cvJ  5cv 

61  =  “  T  +  T  2 

(39) 

6q  —  cv“  —  a6 

(40) 

Fig.  1.  Range  Variance  Versus  F  with  T  —  Is  and  cr„  = 
10  m 


0123456789  IO 

Tracking  Index  r 


Fig.  2.  Correlation  Coefficient  for  Range  and  Range  Rate 
Versus  F  with  T  =  Is  and  crv  ~  10  m 


If  there  is  no  range-Doppler  coupling  (i.e.,  —  =  0)  [3], 
(26)  through  (29),  and  of  (35)  give 


Pk\k  -  cr't 


cv 

0 


T 


r2  = 


0_ 

02(2 1-0) 

2(1  -  a)T2  . 

JL 

1  —  a 


0=  2(2  -  a)  -  4\/l  -  a 


(41) 

(42) 

(43) 


These  results  of  (41)  through  (43)  agree  with  that  derived 
in  the  Appendix  of  [1]  and  Kalata’s  paper  [3]. 

Inserting  (26)  into  (30)  to  eliminate  a  and  expressing  a 
polynomial  equation  of  0  in  terms  of  f,  A<  and  T  gives 


04-l-03  + 


T4 

4 


02  -  iA0  +  r4  =  o 


(44) 


Note  that  (44)  does  not  depend  on  the  sign  of  At.  Newton’s 
method  with  a  good  initial  guess  of  0  can  be  used  to  find  a 
zero  of  (44).  A  good  initial  guess  of  0  is  given  by  the  case 


of  At  =  0,  where  the  initial  0  is  computed  with  (43)  and 
a  is  given  by  [4] 


a  =  - 1  (r2  +  sr  -  (r  +  4)  +  sr)  (45) 

After  0  is  found  by  utilizing  Newton’s  method,  the  steady- 
state  gain,  cv  and  the  error  covariance  are  then  calculated 
At 

for  the  given  T  and  —  using  (26)  through  (29). 

Various  results  for  the  a,0  filter  are  shown  in  Figures 
1  through  7  and  the  results  were  comfirmed  by  iterating 
the  Kalman  filtering  algorithm  until  the  steady-state  con¬ 
ditions  were  achieved.  Figure  1  shows  that  the  variance  of 
the  range  estimate  decreases  with  increasing  At  for  a  given 
F.  Figure  2  shows  that  the  correlation  coefficients  for  range 
and  range  rate  estimates  are  between  zero  and  one.  Note 
that  as  At  becomes  more  positive,  the  correlation  coef¬ 
ficient  for  range  and  range  rate  goes  to  zero  at  higher  F. 
This  means  that  the  range  and  range  rate  estimates  are  less 
correlated  to  each  other  at  higher  F.  Thus  the  accuracy 
of  the  predicted  measurements  will  be  improved  at  more 
positive  At.  Figure  3  shows  that  the  variance  of  range  rate 
decreases  with  increasing  At.  Figure  4  shows  that  cv  de¬ 
creases  with  increasing  At.  For  the  case  of  A t  =  0,  cv  will 
go  to  one  as  F  approaches  infinity.  Figure  5  shows  that  0 


Fig.  3.  Range  Rate  Variance  Versus  V  with  T  =  i  s  and 
c tv  =  10  in 


Fig.  4.  a  Versus  F  with  T  =  Is  and  <jv  =  10  m 

is  maximized  for  no  range- Doppler  coupling  (i.e.,  Ai  =  0) 
and  0  does  not  depend  on  the  sign  of  A t  as  indicated  by 
(44).  For  the  case  of  At  =  0,  0  will  go  to  two  as  F  ap¬ 
proaches  infinity.  In  Figure  6,  0  is  plotted  versus  a  for 
.  At 

various  values  of  Note  that  there  may  be  two  solu¬ 

tions  of  0  for  a  given  a  when  At  >  0,  and  there  is  only  one 
solution  of  0  for  a  given  t*  when  At  <  0.  Figure  7  shows  the 

At 

contours  of  constant  F  for  0  versus  a  for  —1.0  <  —  <1.0. 

At 

Comparing  Figures  6  and  7,  the  curve  for  —  =  0  passes 

through  the  peak  of  the  contours  of  constant  F.  Note  that 
At 

for  —  —  0,  the  values  of  a  and  0  will  go  to  one  and  two 
as  F  approaches  infinity,  respectively. 


IV.  INITIALIZATION  GAINS 


Fig.  6.  0  Versus  a  with  T  =  Is,  crv  —  10  m  and 

At 

0.01  <  F  <  10.0  for  various  values  of  — 


Figure  7.  0  Versus  a  with  T  =  1  s,  crv  =  10  m  and 
At 

— 1.0  <  —  <1.0  for  various  values  of  F 


Least-squares  estimation  can  be  used  to  derive  the 
gains  for  initializing  an  a,0  filter.  For  linear  least-squares 
estimation,  an  equation  formulating  the  measurement  vec¬ 
tor  Z  as  a  linear  function  of  the  parameter  vector  X  to  be 
estimated  is  given  by 

Z  =  WX  -f  V  (46) 

where  E[V]  =  0  and  E[V  V1]  =  cr'^Ijsf  for  N  +  1  measure¬ 
ments  with  lN  denoting  the  (N  +  1)  x  (N  +  1)  identity 


matrix.  The  least-squares  estimate  X  [4]  is  given 

X  =  (WTW)~1WTZ  (47) 

with  error  covariance 

PN  =  COV[X]  =  <tI{WtW)~1  (48) 

A  current  state  estimate  of  X  denoted  by  Aq  — 
[r0  f[)]r  can  be  obtained  from  the  current  measurement 


plus  the  N  previous  measurements  as 

Z\i  =  U'-vAo  +  Vn  (49) 


where 


Z,v  =  [y-;V  V-N+ 1  ^-1  yu] 

V'v  =  [tr-iV  +  l  i;-l  Vo] 

1  1  ...  1 
At -NT  At  +  (— iV  +  1)T  ...  At 


(50) 

(51) 

(52) 


with  A'o|-i  =  [0  Op  .  Note  that  as  k  approaches  infinity, 
the  functions  of  k  for  ak  and  0k  go  to  zero.  This  is 
because  the  least-squares  estimation  assumes  there  is  no 
process  noise.  However,  the  steady-state  gains,  a  and 
0  are  nonzero  in  the  Kalman  filter,  which  takes  process 
noise  into  consideration.  If  aLV  =  0  (e.g.,  exoatmospheric 
ballistic  target),  and  a  =  0  —  0,  c\k  and  0k  are  valid  for 
longer  periods.  Using  the  range  measurement  at  t  =  0,  yo, 
and  the  range  measurement  at  t  —  T,  \j\ ,  the  initial  state 
estimate  and  error  covariance  based  on  two  measurements 
with  an  LFM  waveform  are  given  by 


where  £*[lqv]  =  0.  E'[V.\V$]  =  (t-J/jv,  At  is  the  range- 
Doppler  coupling  coefficient  in  (14)  and  Wqv  is  a  combina¬ 
tion  of  the  state  transition  matrix  and  the  output  matrix. 
The  error  covariance  of  the  least-squares  state  estimate  Xy 
with  an  LFM  waveform  is  given  by 

CW[.Y„]  -  Cv  =  t%(WlWN)~x  (53) 

where 

At  _  N  T 
;V  +  1  N  + 1  2 

( Afc)“  — At;VT  .  N(‘2N 4-1 )  TJ 
/V+l  +  N ■+*  1  6  J 

(54) 


The  A'i|i  and  Ppi  are  the  result  of  the  least-squares 
estimation  with  the  first  two  measurements. 


V.  CONCLUDING  REMARKS 


(WlWN) 


At  - 

N  +  l  2 


Then 


Tn  = 


2(ri! 


2;V  +  1 


At  6 

bT+N 


T 


(1 


2  At 
NT 


Ti-— ) 

Tv  ;VT' 


6 


NT2 


(55) 


where  c//v  =  (iV  +  l)(iV  +  2).  Note  that  At  cannot  be 
selected  to  minimize  the  determinant  of  |/V|,  because 

At  N 

|/^|  is  not  a  function  of  At.  Also  note  that  for  —  =  —  , 

the  variance  of  range  is  minimized  and  the  range  and  range 
rate  estimates  are  uncorrelated. 

Let  the  scheduled  gains  for  initialization  for  the  a,  0  filter 
be  denoted  by 


K„  =  [  ak  T]T  (56) 

where  k  denotes  the  gains  for  processing  the  k  4-  Ist  mea¬ 
surement. 

Letting  Pk\k  —  I\  of  (55)  in  (8)  provides  a  simple  gain 
scheduling  procedure  for  a  and  0  during  initialization  that 
is  given  by 


2(2fc  +  l)-6—  . 

afc  =  rnan  {k  +  1) (Ar  +  2)  ’ 

r  6  ^ 

A= '"“{(*  +  i)n  +  2) -0) 


(57) 

(58) 


The  steady-state  gains  and  error  covariance  for  the  a. 
0  filter  are  presented  as  a  simple  function  of  a,  0 ,  At ,  T  and 
a2.  The  results  for  the  cv,  0  filter  were  confirmed  with  the 
Kalman  filtering  equations.  The  initialization  equations 
of  (60)  and  (61)  assume  two  identical  LFM  waveforms 
at  t  =  0  and  t  =  T.  However,  in  a  phased  array  radar, 
the  confirmation  dwell  may  include  a  LFM  waveform  with 
up-chirp  and  down-chirp  to  achieve  range  and  range  rate 
measurements  from  a  single  dwell. 
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Steady-State  Tracking  with  LFM  Waveforms 


Correspondence 

The  steady-state  gains  and  error  covariance  are  derived  for 
a  two-state  Kalman  filter  (i.e.,  an  a,  /3  filter)  for  tracking  with 
linear  frequency  modulated  (LFM)  waveforms.  A  procedure  is 
given  for  calculating  the  a,  j3  gains  from  the  tracking  index  F , 
the  sample  period  7,  and  the  range-Doppler  coupling  coefficient 
Ar.  The  steady-state  error  covariance  is  found  to  be  a  simple 
function  of  a,  0,  At,  T  and  the  measurement  noise  variance 
The  expressions  for  the  steady-state  gains  and  error  covariance 
were  confirmed  numerically  with  the  Kalman  filtering  equations. 
A  gain  scheduling  technique  for  a  and  f3  during  initialization  is 
also  given. 


INTRODUCTION 

The  Kalman  filter  produces  an  optimal  estimate 
of  the  target  state  when  given  the  motion  model  and 
a  sequence  of  sensor  measurements  corrupted  with 
white  Gaussian  errors.  The  computational  burden  of 
maintaining  the  Kalman  filter  may  prohibit  its  use 
when  many  targets  are  being  tracked.  Steady-state 
filters  are  also  good  for  filter  design  and  analysis. 
Constant  data  rate  and  constant  covariances  required 
for  steady-state  filters  are  not  typically  true  in  real 
systems,  so  approximate  gains  are  typically  used. 

The  gains  are  based  on  the  steady-state  gains  of 
the  Kalman  filter  and  implemented  either  using  a 
fixed  gain  or  an  easily  generated  gain  schedule. 

The  resultant  filter  for  tracking  range  and  range 
rate  is  called  the  a,  0  filter,  where  approximate 
filter  gains  a,  0  may  be  used  instead  of  the  optimal 
Kalman  gains  [1].  In  [2],  the  range-Doppler  coupling 
of  linear  frequency  modulated  (LFM)  waveforms 
is  shown  to  have  a  significant  effect  on  tracking 
accuracy.  In  [3]  and  the  Appendix  of  [1],  the 
steady-state  gains  and  error  covariance  are  derived 
for  range  measurements  without  range-Doppler 
coupling.  Here,  the  range-Doppler  coupling 
associated  with  LFM  waveforms  is  included  in 
the  measurement  equation,  and  expressions  for 
the  steady-state  gains  and  error  covariance  are 
calculated. 
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MATHEMATICAL  MODEL 

The  dynamics  model  commonly  assumed  for  a 
target  in  track  is  given  by 

Xk+i  =FkXk  +  Gkwk  O) 

where  wk~N(0,Qk)  is  the  process  noise,  Fk  defines 
a  linear  constraint  on  the  dynamics,  and  Gk  is  the 
input  matrix  for  maneuver.  The  target  state  vector  Xk 
contains  the  range  and  range  rate  of  the  target  at  time 
k.  The  linear  measurement  model  is  given  by 

Yk  =  HkXk  +  vk  (2) 

where  Yk  is  typically  the  target  position  measurement 
and  v*  ~  N(0,Rk)  is  the  measurement  noise,  and  Hk 
is  the  output  matrix  that  includes  range-Doppler 
coupling  coefficient. 

The  measured  range  and  the  variance  for  an  LFM 
waveform  and  the  range-Doppler  coupling  coefficient 
for  an  LFM  waveform  are  given  by  [2]  as 


yk  =  rk  +  Atrk  +  vk 

(3) 

2  nr  2l  A t2c2  c 2 

<rv-E[vk\  -  g£r2 -  8B2sr 

(4) 

A,  =  /»"  =  *T 

f2~fl  & 

(5) 

where  rk  is  the  true  range  at  time  k ,  5ft  is  the 
signal-to-noise  ratio,  c  is  the  speed  of  light,  f0 
is  the  nominal  carrier  frequency,  fx  is  the  initial 
frequency,  and  f2  is  the  final  frequency.  Also  \B\  is  the 
bandwidth  of  the  LFM  waveform  and  r  is  the  pulse 
length.  Note  that  the  demodulation  frequency  can  be 
modified  with  the  range  rate  estimate  to  reduce  the 
bias  in  the  range  measurement  due  to  range-Doppler 
coupling.  However,  if  the  range  rate  estimate  is 
used  to  remove  the  bias  in  the  range  measurements, 
the  measurement  errors  become  correlated  with 
the  state  estimate  and  the  typical  Kalman  filter 
must  be  modified  to  account  for  the  state  correlated 
measurement  errors. 

The  Kalman  filtering  equations  associated  with  the 
state  model  in  (1)  and  the  measurement  model  in  (2) 
are  given  by  the  following  equations. 

Time  Update : 

%k\k- 1  =  ^-l^-l|fe-l  W 

Pk\k-\  ~  Fk-\Pk-\\k-\Fk-\  +  ^k-\Qk-\^k-\ *  0) 

Measurement  Update : 

Kk  =  +  Rk1  (8) 

=  +  Xk\?k  ~  HkXk\k-l]  @) 

Pk\k  =  U-KkHk)Pk]k_l  (10) 

where  the  subscript  notation  ( k  \  j )  denotes  the  state 
estimate  for  time  tk  when  given  measurements  through 


time  tj.  The  Kk  is  the  Kalman  gain  that  minimizes 
the  mean  square  error  in  the  state  estimate.  Also, 
note  that  Xk  ~  N(Xk\k,Pk]k)  with  XkVc  and  Pk\k  denoting 
the  mean  and  error  covariance  of  the  state  estimate, 
respectively. 


THE  a,P  FILTER 


For  the  Kalman  filter  in  steady-state  conditions, 

Pk\k  =  3-iit-i.  md  pk+i\k  =  30(1  Kk  =  Kk- 1;  For  a 

Kalman  filter  to  achieve  these  steady-state  conditions, 
the  error  processes,  wk  and  v^,  must  have  stationary 
statistics  and  the  data  rate  must  be  constant.  When 
the  noise  processes  are  not  stationary  or  the  data 
rate  is  not  constant,  a  filter  using  the  steady-state 
gains  provides  suboptimal  estimates.  The  a,P  filter 
is  the  steady-state  Kalman  filter  for  tracking  nearly 
constant  range  rate  targets.  The  a,  P  filter  .is  a  single 
coordinate  filter  that  is  based  on  the  assumption 
that  the  target  is  moving  with  constant  range  rate 
plus  zero-mean,  white  Gaussian  acceleration  errors. 
Given  this  assumption,  the  filter  gains  a  and  p  are 
chosen  as  the  steady-state  Kalman  gains  that  minimize 
the  mean-square  error  in  the  range  and  range  rate 
estimates. 

For  the  a,P  filter  for  piecewise  constant 
acceleration  errors, 


II 

(ii) 

■l  ri 

(12) 

Fk  =  F  = 

.0  lj 

II 

li 

o 

T2  ]T 

h  T\ 

(13) 

Hk  =  H  =  (1  Af] 

(14) 

Rk  =  R  =  <tv2 

(15) 

Qk  =  Q  = 

(16) 

II 

II 

p]T 

la  rJ 

(17) 

where  rk  and  rk  are  the  range  and  range  rate  of  the 
target,  respectively,  a2  and  c2  are  the  variances  of 
measurement  noise  and  process  noise,  respectively,  T 
is  the  sample  period  between  measurements,  At  is  the 
range-Doppler  coupling  coefficient,  and  a  and  P/T 
are  the  steady-state  Kalman  filter  gains. 

Let  the  steady-state  error  covariance  matrix  of  the 
filtered  estimates  for  the  a,  p  filter  be  denoted  as 


^t+i|*+i  -  Pk\k~  P  - 


P  n 
■Pl2 


P 12 
P22k 


(18) 
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Using  (14),  (15),  and  (18)  in  (8)  gives  the  steady-state 
gain  as 


Kk  =  Kk_l  =K  =  PHJR~l  = 


T  n- 1  [Pu+Pn&t}  -1 

R  =  av  . 

■P 12  +  Pn^- 


Using  (17)  with  (19)  gives 

a  =  On  +PnAt)av2  ( 

P  =  (Pl2  +  P22  At)T(T~2.  ( 

Inserting  (7)  into  (10)  and  setting  Pk\k  = 

P  for  steady-state  conditions  gives 


P<r$(2a-p) 
2T2  fl  —  a- 


Substituting  pl2  into  (20)  gives 


Pu  —  aav 


1  — a  — /?- 


x  1  -a- 


2 a  +  P\  Arl  At 


T  T 


[I  -  KH ]  lP  =  FPFr  +  GGra2 .  (22)  The  steady-state  error  covariance  is  then  given  by 


Pk\k  -  av 


i  rtAr 

1  —  a  —  0— 
T 


(  At 

r  f 1 —  a  —  0— 


2 a  +  p\  Arl  At 

2~~  JT  T 


2a  +  p  \  Ar 

2  J  7" 


rfl-a-/3y 


1  -a  — 


2d  +  P  \  A? 


P(2a-p) 

T*(l-a-0j- 


•  (30) 


[/-Z/f]"1?  = 


FPF1  = 


pn+aAtpn 


jPn  +  (1  a)Pi2 


PH  +  2Tpn  +  r2p22  Pl2  +  TP22 
Pl2  +  TP22  P22 

r  r4  r3i 


cc^i-oi  2  . 

Li-  H 


,  PAt\  . 

I  —  J  /?12  +  a  Ar/?22 


"P12  +  (1  ~  a)P22 


Equating  the  (1,2)  and  (2,2)  elements  of  (22)  gives  Equating  the  (1,1)  elements  of  (22)  and  using  (27) 

expressions  for  p12  and  p22.  Substituting  p12  and  p22  through  (29)  gives 

into  (21)  gives  .  ~ 

2  2  a3/3i+a2p2  +  a1/3  +  a0  =  0  ( 

t^2  _  1  aw  _  P  ,*>c\ 


(26)  where 


1  -0-/3- 


where  T  is  the  tracking  or  maneuvering  index.  Using 
(26)  to  eliminate  <7^,  in  the  expressions  of  pl2  and  p22 
gives 


T  1  -a-p- 


2 a  +  p\  At 


A  t\  A  t 


/ArV  Ar  Ar  1-a 
a2  =  “2a  I  —  1  -2ay+3y  +  -^— 

=  —3a2  ^  —  a2  +  2a^  +  3a  —  2 


<zn  =  a  —  a  . 
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There  are  three  solutions  for  (31).  By  comparing  with 
the  Kalman  filter  numerical  results,  the  valid  solutions 
of  (31)  for  P  are  given  by 


(36) 


For  At/T  >  0,  both  solutions  of  (36)  are  valid  for 
p  depending  on  the  selection  of  F  and  At/T.  The  p_ 
(i.e.,  minus  sign  in  the  front  of  the  radical)  is  used  for 
P  with  increasing  a  until  the  maximum  a  is  reached, 
which  is  given  by 


Then  P+  (i.e.,  plus  sign  in  the  front  of  the  radical)  is 
used  for  P  from  amax  with  decreasing  a  and  it  is  valid 
as  long  as  1  -  a  -  P+At/T  >  0,  that  is,  F  in  (26)  is  a 
real  number.  For  At/T  <  0,  3__  is  the  only  solution  for 
p.  If  At/T  =  —  f,  (31)  becomes  a  quadratic  equation 
given  by 

b2p2  +  blp  +  b0  =  0  (38) 


where 


h  -  a  1 

i>T2 


b>  =  -T  +  T-2  <40 

b0  =  a2-a3.  (41 

If  there  is  no  range-Doppler  coupling  (i.e.,  At/T  =  0) 
[1,  3,  7],  (30),  (26)  and  p_  of  (36)  give 

r  @  ~\ 

_  ,f“  r 


Pk\k-°*  p  p2(2a-p)  ( 

-T  2(1 -a)T2- 

r2  =  -  JL.  (43 

o2  1  -  a 

P  =  2(2  —  a)  —  A\J  1  —  a.  (44 

These  results  of  (42)-(44)  agree  with  that  derived 
from  [1,  Appendix]  and  Kalata’s  papers  [3,  4]. 

Inserting  (26)  into  (31)  to  eliminate  a  and 
expressing  a  polynomial  equation  of  P  in  terms  of  F , 
At  and  T  gives 


p*-r2p2  + 


■  2T2  F02-rV  +  r4  =  o. 


Note  that  (45)  does  not  depend  on  the  sign  of  At. 
Newton’s  method  with  a  good  initial  guess  of  P  can 
be  used  to  find  a  zero  of  (45)  when  F  and  At/T  are 
given.  A  good  initial  guess  of  P  is  given  by  the  case 
of  At  =  0,  where  the  initial  P  is  computed  with  (44) 
and  a  is  given  by  [6] 

a  =  -i(r2  +  sr  -  (r  +  4)Vr2  +  8r).  (46) 

After  P  is  found  by  utilizing  Newton’s  method,  the 
steady-state  gain,  a,  and  the  error  covariance  are  then 
calculated  for  the  given  T  and  At/T  using  (26)  and 
(30).  These  results  were  confirmed  with  the  Kalman 
filtering  equations. 

Various  results  for  the  a,P  filter  are  shown 
in  Figs.  1-7  and  the  results  were  comfirmed  by 
iterating  the  Kalman  filtering  algorithm  until  the 
steady-state  conditions  were  achieved.  Fig.  1  shows 
that  the  variance  of  the  range  estimate  decreases 
with  increasing  At  for  a  given  T .  Fig.  2  shows  that 
the  correlation  coefficients  for  range  and  range  rate 
estimates  are  between  zero  and  one.  Note  that  as  At 
becomes  more  positive,  the  correlation  coefficient  for 
range  and  range  rate  goes  toward  zero  at  higher  F . 

This  means  that  the  range  and  range  rate  estimates 
are  less  correlated  with  each  other  at  higher  F .  Thus 
the  accuracy  of  the  predicted  measurements  will 
be  improved  at  more  positive  At.  Fig.  3  shows  that 
the  variance  of  range  rate  decreases  with  increasing 
At.  Fig.  4  shows  that  a  decreases  with  increasing 
At.  For  the  case  of  At  =  0,  a  will  go  to  one  as  T 
approaches  infinity.  Fig.  5  shows  that  3  is  maximized 
for  no  range-Doppler  coupling  (i.e..  At  =  0)  and  P 
does  not  depend  on  the  sign  of  At  as  indicated  by 
(45).  For  the  case  of  At  =  0,  P  will  go  to  two  as  F 
approaches  infinity.  In  Fig.  6,  p  is  plotted  versus  a  for 
various  values  of  At/T.  Note  that  there  may  be  two 
solutions  of  p  for  a  given  a  when  At  >  0,  and  there 
is  only  one  solution  of  P  for  a  given  a  when  At  <  0. 
Fig.  7  shows  the  contours  of  constant  F  for  p  versus 
a  for  -1.0  <  At/T  <  1.0.  Comparing  Figs.  6  and  7, 
the  curve  for  At/T  =  0  passes  through  the  peak  of 
the  contours  of  constant  F .  Note  that  for  At/T  =  0, 
the  values  of  a  and  P  will  go  to  one  and  two  as  T 
approaches  infinity,  respectively. 

INITIALIZATION  GAINS 

Least-squares  estimation  can  be  used  to  derive 
the  gains  for  initializing  an  a,P  filter.  For  linear 
least-squares  estimation,  an  equation  formulating 
the  measurement  vector  Z  as  a  linear  function  of  the 
parameter  vector  X  to  be  estimated  is  given  by 

Z  =  WX  +  V  (47) 

where  £[V]  =  0  and  E[VVT ]  =  a2lN  for  N  +  1 
measurements  with  IN  denoting  the  ( N  +  1)  x  (V  +  1) 
identity  matrix.  The  least-squares  estimate  X  [5]  is 
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Tracking  Index  r 


Fig.  1.  Range  variance  versus  T  with  T  =  1  s  and  av  =  10  m. 


Fig.  2.  Correlation  coefficient  for  range  and  range  rate  versus  F  with  T  =  Is  and  <rv  =  10  m. 


given 

measurement  plus  the  N  previous  measurements  as 

X  =  (WTW)~XWTZ 

(48) 

ZN=WNX0  +  VN 

I 

with  error  covariance 

PN  =  cov[X]  =  a2v(WTWrl. 

A  current  state  estimate  of  X  denoted  by 

X0  =  [r0  r0]r  can  be  obtained  from  the  current 

(49) 

where 

ZN  =  [y_N  y_^+1  •*'  i 

=  V_N  +  1  •••  V_! 

J'oi  1 

I'q]  1 

CORRESPONDENCE 


Tracking  Index  r 

Fig.  3.  Range  rate  variance  versus  T  with  T  =  Is  and  crv  =  10  m. 


Tracking  Index  r 

Fig.  4.  a  versus  T  with  T  =  Is  and  <rv  =  10  m. 


A  t-NT  Al  +  (-N+l)T  ■■■  At-T  At 


(53) 

where  E[VN]  =  0,  =  a2vlN,  At  is  the 

range-Doppler  coupling  coefficient  in  (14)  and 
WN  is  a  combination  of  the  state  transition  matrix 
and  the  output  matrix.  The  error  covariance  of  the 
least-squares  state  estimate  X0  with  an  LFM  waveform 


is  given  by 

cov[X0]  =  PN  =  o2v(W*WNrl  (54) 

where 

'  1  M-N^ 

(WT\Vn)  =  (N+1) 

.At-N-  (At)2  -  AtNT  +  N(2N  +  1)— 

2  d 

(55) 
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Tracking  Index  T 


Fig.  5.  p  versus  T  with  T  =  Is  and  crv  =  10  m. 


Fig.  6.  p  versus  a  with  T  =  1  s,  av  =  10  m,  0.001  <  P  <  10.0  for  various  values  of  A t/T. 


Then 


2cr2 


cov[X0]  PN-  (Ar+1)(^  +  2) 


At 

2N  +  1  -6-jT  + 


6^/Ar\2  3  /  2Af> 

N  \  T  )  TV  NT  j 


of  At.  Also  note  that  for  A  t/T  =  N/ 2,  the  variance 
of  range  is  minimized  and  the  range  and  range  rate 
estimates  are  uncorrelated. 

Let  the  scheduled  gains  for  initialization  for  the 
a,/3  filter  be  denoted  by 


3  ( 

2AA 

6 

1 

J 

T  \ 

NT) 

NT2  J 

(56) 

ii 

C 


Note  that  At  cannot  be  selected  to  minimize  the 
determinant  of  PN,  |/^|,  because  \PN\  is  not  a  function 


where  k  denotes  the  gains  for  processing  the  k  +  1st 
measurement. 


CORRESPONDENCE 


Letting  tfc\k  ~  of  (56)  in  (8)  provides  a  simple 
gain  scheduling  procedure  for  a  and  /?  during 
initialization  that  is  given  by 


olk  ~  max  < 


2(2* +  1)- 6^ 

(*  +  l)(*  +  2) 

6 


ft  =  nm{(it7WT2)’g} 


(58) 

(59) 


with  X0|_j  =  [0  Of.  Note  that  as  k  approaches 
infinity,  the  functions  of  k  for  ak  and  (3k  go  to  zero. 
This  is  because  the  least-squares  estimation  assumes 
there  is  no  process  noise.  However,  the  steady-state 
gains,  a  and  /3,  are  non-zero  in  the  Kalman  filter, 
which  takes  process  noise  into  consideration.  If  aw  =  0 
(e.g.,  exoatmospheric  ballistic  target),  and  a  =  (3  = 

0,  ak  and  pk  are  valid  for  longer  periods.  Using 
the  range  measurement  at  t  =  0,  y0,  and  the  range 
measurement  at  t  =  7\  y{7  the  initial  state  estimate  and 
error  covariance  based  on  two  measurements  with  an 
LFM  waveform  are  given  by 


The  X1(1  and  P1(1  are  the  result  of  the  least-squares 
estimation  with  the  first  two  measurements. 

CONCLUDING  REMARKS 

The  steady-state  gains  and  error  covariance  for 
the  a,  (3  filter  are  presented  as  a  simple  function  of 
a,  /?,  At,  T,  and  cr^.  The  results  for  the  a,/3  filter 
were  confirmed  with  the  Kalman  filtering  equations. 
The  initialization  equations  of  (60)  and  (61)  assume 
two  identical  LFM  waveforms  at  t  =  0  and  t  =  T. 
However,  in  a  phased-array  radar,  the  confirmation 
dwell  may  include  an  LFM  waveform  with  up-chirp 
and  down-chirp  to  achieve  range  and  range  rate 
measurements  from  a  single  dwell. 
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Introduction 


•  Fusion  is  only  a  supporting  process  and  it  is  unlikely  that  any 
customer  is  going  for  pay  for  “fusion  for  the  sake  of  fusion.” 

•  In  many  typical  systems,  the  developers  of  the  components  (e.g., 
sensors)  and  system  integrators  most  often  do  not  understand 
fusion  algorithms  nor  appreciate  their  value. 

•  Critical  to  develop  and  apply  fusion  algorithms  in  a  systems 
context  and  optimize  the  system-level  performance. 

•  Demonstrations  of  the  potential  for  improvements  in  system-level 
performance  will  encourage  system  developers  to  take  the  risk 
associated  with  implementing  fusion  algorithms. 

•  Promising  areas  for  potential  improvements  in  system-level 
performance  of  sensor  systems  is  the  use  fusion  for  efficient 
utilization  of  sensor  resources  and  control  of  sensor  modes. 
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Introduction 


•  Tracking  benchmark  is  a  set  of  scenarios  and 
associated  simulation  environment  that 

>  Captures  salient  features  of  the  tracking  system  of  interest 

>  Provides  a  “level  playing  field”  for  evaluation  and  comparison 
of  candidate  tracking  algorithms 

•  Expected  benefits  of  benchmark  process 

>  Promote  research  and  development  with  respect  to  many 
“real-world”  issues  associated  with  the  deployment  of  the 
tracking  system  of  interest 

>  Development  of  algorithms  that  optimize  system-level 
performance  metrics  of  the  tracking  system  of  interest 

Georgia  J,  [stesaPsOn 

- - - - -  Tech  U  OmsSSKMis 


Tracking  Requirements 


Tracking  In  Legacy  Simulations 


NSWC  Tracking  Benchmark  I 


Average  Sample  Periods  versus  Computational  Cost 
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Performance  Metrics 


•  Completeness  history 

-  Composite  completeness  (implemented) 

-  Relative  completeness 

•  Timeliness 

-  Composite  track  initiation  time  (implemented) 

-  Composite  track  relative  initiation  time 

-  Composite  track  convergence  time 

•  Track  continuity 

-  Cumulative  swaps  of  composite  tracks  (implemented) 

-  Cumulative  broken  composite  tracks  (implemented) 

•  Ambiguity 

-  Composite  track  redundant  track  mean  ratio  (implemented) 

-  Composite  track  spurious  track  mean  ratio  (implemented) 

•  Accuracy 

-  Composite  track  accuracy  (implemented) 

-  Composite  track  covariance  consistency  (implemented) 

•  Cross-platform  commonality  history 

-  Ratio  of  non-common  composite  track  numbers  (implemented) 

-  Composite  track  state  estimate  differences 

•  JCTN  loading 

-  Communication  data  loading  (implemented)  i 

-  Processor  loading  Georgia  i  Dtowmafo 
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Reporting  Responsibility 
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Concluding  Remarks 


Three  areas  critical  to  sensor  fusion  are: 


•  Optimization  of  system-level  performance  through  the  use 
of  fusion 

•  Extension  of  fusion  algorithm  to  address  sensor 
phenomenology 

•  Characterization  of  the  performance  of  different  fusion 
algorithms  versus  the  cost  of  implementation  via 
benchmark  problems 
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Introduction 


Multiplatform-Multisensor  Architectures 
Coordinate  Systems 

Multiplatform-Multisensor  Track  Filtering 

-  Target  Motion  Modeling 

-  Radar  Measurements 

-  Extended  Kalman  Filter 

-  Processing  Local  Measurements 

-  Processing  Remote  Measurements 

-  Processing  Time-delay,  Out-of-Sequence  Measurements 

-  Interacting  Multiple  Model  (IMM)  Estimator 


Conclusion  Remarks 
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^ntmductiorj 

Target  Tracking  includes  two  types  of  uncertainty 

-  Target  Motion  Uncertainty:  magnitude  and  mode 

-  Measurement  Uncertainty:  magnitude  and  mode 

Local  Measurement  refers  to  any  sensor  measurements  that  originate  on  the  platform 
where  the  track  processor  resides. 

Remote  Measurement  refers  to  any  sensor  measurements  that  originate  on  a  platform 
other  than  the  platform  where  the  track  processor  resides. 

Local  Track  refers  to  any  track  that  includes  only  local  measurements. 

Remote  Track  refers  to  any  track  that  includes  only  remote  measurements. 

Composite  Track  refers  to  any  track  that  includes  measurements  from  more  than  one 
sensor. 


Georgia  Tech  Research  Institute 

Introduction 

Gridlock  refers  to  the  process  by  which  data  from  one  frame  of  reference  is  transformed  to 
another  frame  of  reference.  The  frame  of  reference  may  be  the  tracking  frame  of  reference 
and/or  the  sensor  frame  of  reference.  This  process  will  include  a  geometric  transformation 
component  and  an  intersensor  bias  estimation  and  removal  component.  The  term  gridlock 
may  be  used  interchangeably  with  sensor  data  registration. 

Platform-centric  Tracking  refers  to  those  tracking  systems  where  the  sensors  are  all 
located  at  a  single  site  or  on  a  single  platform  (c.g..  a  ship  or  an  aircraft).  In  these  cases, 
the  transformation  of  data  from  one  frame  of  reference  to  another,  which  is  the  first  step  in 
the  gridlock  process,  can  be  accomplished  by  rather  simple  transformations  that  change 
very  slowly  with  time.  The  data  passing  from  the  various  sensors  to  the  tracker  can  he 
accomplished  in  a  relatively  simple  manner  as  well  (e.g.,  over  a  local  area  network). 

Network-centric  Tracking  refers  to  those  tracking  systems  where  sensors  may  be 
located  at  significant  distances  from  one  another  as  in  the  case  of  radar  data  from  multiple 
ships  or  a  ship  and  an  aircraft.  In  these  cases,  the  transformation  of  data  from  one  site  to 
another  is  more  involved,  and  the  existence  of  a  communication  link  is  implied. 
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Figure  1:  Illustration  of  network-centric  tracking 
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Figure  2:  Assigned  reporting  responsibility 
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Architectures! 


Table  1:  Assumptions  for  Hypothetical  Network- Centric  Tracking  Systems 


Max  position  =  500  km 

Position  LSB  =  lm 

Max  velocity  =  1,000  m/s 

Velocity  LSB  =  0.1  m/s 

Max  range  =  500  km 

Range  LSB  =  lm 

Range  accuracy  =  1  -  300m 

Max  bearing  —  360° 

Bearing  LSB  =  0.05° 

Max  elevation  =  90° 

Elevation  LSB  =  0.05° 

Bearing  accuracy  =  0.05  -  5° 

Elevation  accuracy  =  0.05-5° 

Max  time  —  86,400  s 

Time  LSB  =  0.1  s 

Max  track  number  (ATtrac ks)=  1-000 

Track  number  LSB  =  1 

Number  of  platforms  =  5 

Tracking  interval  (/track)  =  5  s 

Georgia  Tech  Research  Institute 


Assigned  Reporting  Responsibility! 

Number  of  bits  required  to  represent  each  quantity  can  be  computed  as  follows: 


Velocity  LSB 
Max  time 


J  =>  ?7-vei  =  14  bits  4*  1  sign  bit 


Time  LSB 


)  =>  ntimc  =  20  bits 


(1) 

(2) 

(3) 

(4) 


Thus,  20  bits  are  necessary  for  each  position  quantity,  15  bits  for  each  rate  (velocity)  quantity, 
20  bits  for  the  time  quantity,  and  10  bits  for  the  track  number.  Thus,  the  total  data  rate 
required  to  support  the  hypothetical  system  is  given  by 


^pos 

> 

log2( 

^vel 

> 

log2( 

W'time 

> 

l0§2( 

^number 

> 

log2( 

/foistRR  — 


A^t  racks 
It  rack 


(3npos  -f-  377.Vel  4"  Ti,  time  4“  ^number) 


(5) 


which  results  in  a  data  rate  requirement  of  27  kbps. 
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Figure  3:  Central-level  composite  tracking:  pure 
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Figure  4:  Central-level  composite  tracking:  practical 
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Central-level  Composite  Tracking:  Practical^ 


Data  rate  required  to  pass  just  the  measurements  and  associated  accuracies  with  a  valid  time 
tag  is  computed  with 


Grange 

—  19  bits 

flr_acc  —  ^  bits 

(6) 

^bearing 

=  13  bits 

ttb-acc  =  7  bits 

(7) 

^-elevation 

—  12  bits 

(including  a  sign  bit) 

(8) 

^e_acc 

—  7  bits 

ntime  =  20  bits 

(9) 

For  this  approach,  the  data  rate  required  to  pass  all  measurements  from  legitimate  target 
returns  is  given  by 


^mcas 

^mcas_acc 

-^practical  central 


Grange  “I"  ^bearing  T*  ^elevation 

^■r_acc  d~  n,k_acc  +  n0_acc 
^platforms  ^meas /scan 


^track 


(^mcas  T  ^moas_acc  T  rifime) 


(10) 

(11) 


(12) 


To  support  the  passing  of  AMRs  associated  with  the  maximum  number  of  targets 
(Nmcas/scan  =  1,000),  the  data  rate  requirement  is  17.4  kbps  per  platform  and  87  kbps  for  the 
entire  five-platform  network. 
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Figure  5:  Distributed  tracking  with  track  fusion  for  composite  tracking:  central-level  track  fusion 
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Central-level  Track  Fusion 


Continuing  the  example  from  above,  with  the  added  assumption  that  the  elements  of  the 
covariance  matrix  will  all  be  represented  as  32-bit  quantities,  the  data  rate  requirement  can  be 
estimated  as 


•Rtrack 


’  platforms  ^meas/scan 
-^track 


X  (3  71 

pos  3nvel  “1“  -^number  H"  ^time  T  21nacc) 


where  nacc  is  the  number  of  bits  required  to  accurately  represent  the  elements  of  the  covariance 
matrix.  The  minimum  data  rate  required  to  pass  the  track  data  is  161.4  kbps  per  platform  or 
807  kbps  for  the  five-platform  network. 
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Figure  7:  Distributed  composite  tracking 
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Distributed  Composite  Tracking^ 


The  data  rate  for  passing  measurement  and  accuracy  data  for  the  distributed  composite 
tracking  approach  is  given  by 


-^dist 


-^platforms  ^meas/scan 


composite  r 

-‘track 

X  (^mcas  ~b  ^meas_acc  "b  Tl number  “b  ^time) 


(14) 


Using  the  assumptions  for  the  previous  examples,  the  data  rate  for  this  approach  is  determined 
to  be  19.4  kbps  per  unit  or  97  kbps  for  the  five-platform  network. 
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Table  2:  Data  Rate  Requirements  for  Hypothetical  Network-Centric  Tracking  Systems 


AMR  or 

Accuracy 

Total  rate 

Total  rate 

Architecture 

state  vector 

Update 

Update 

Unit 

Network 

(bits) 

(bits) 

(kbps) 

(kbps) 

Reporting 

responsibility 

135 

NA 

NA 

27 

Cent,  tracking: 

Pure 

- 

- 

- 

- 

Practical 

64 

23 

17.4 

87 

Dist.  tracking: 

Cent.  Track  Fusion 

135 

672 

161.4 

807 

Dist.  Track  Fusion 

135 

672 

161.4 

807 

Tracklets  (4:1) 

135 

672 

40.4 

202 

Composite  Track 

74 

23 

19.4 

97 
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Coordinate 


System^ 


The  defining  parameters  are  the  length  of  the  semi-major  axis,  a.  of  the  generating  ellipse  and, 
either  the  length  of  the  semi-minor  axis.  b.  the  eccentricity,  e,  where 


b 2 


e  = 


or  the  flattening 


/ 


a  —  b 


(15) 


(16) 


The  ellipsoid  is  a  very  good  fit  to  the  gcoid;  the  root-mcan-square  (RMS)  deviation  of  the  gcoid 
from  the  ellipsoid  is  30.5  m  based  on  a  1°  by  1°  worldwide  grid. 
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Coordinate  Systems 


Table  3:  WGS-84  Parameters 


Parameter 

Value 

semi-major  axis  (a) 

eccentricity  (e) 

flattening  (/) 

Earth’s  rotation  rate  (cje) 

6,378,137m 

0.08181919 

1 

298.2572 

7.292115  x  10~5  rad/s 

^Cood|nate^stemsj| 


Position  vector  Ve  to  a  point  above  Earth’s  surface  can  be  specified  in  frame  Te  as 

{N  +  H)  c(f>  cl 
(N  +  H)  c (f>  s l 
(iV(l-e2)  +  K)  s <P 


(17) 


where  <j>  is  the  geodetic  latitude.  £  is  the  geodetic  longitude,  T~L  is  the  geodetic  height,  N  is  the 
prime  vertical  radius  of  curvature  (PVRC). 

The  PVRC  is  a  function  of  the  geodetic  latitude  and  is  calculated  from 

a 


N  = 


y/l  -  e 2  s 2<j) 


(18) 


where  a  is  the  length  of  the  semi-major  axis  of  the  WGS-84  generating  ellipse  and  e  the 
eccentricity. 
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Coordinate  Systems 


The  composite  rotation  or  direction  cosine  matrix  from  ECEF  to  frame  m  is  given  by 


Cm|e  =  -  s<j)m  ctm  -  s <f>m  s£m  c.(j)n 


C0m  ct 


C (j>m  S im  S (j)v 


where  (j)m  and  tm  arc  the  latitude  and  longitude  of  frame  m. 

The  position  of  target  T  relative  to  platform  n  is  denoted  by  position  vector  lZn .  The 
corresponding  position  vector  from  platform  rn  is  given  by 

7Zm  =  Cm]n1Zn  +  Dm|n 


where 


£im\n  _  £»m|c£»e!n  _  —  1  _  ^-j7n\c^Qn\c^T 

£>m|n  =  Cmle(V'-Vem) 
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Coordinate  Systems  | 

The  elements  of  Cm^n  are  given  by 

Cm|n(  1,1) 

-  cAf”1” 

(23) 

C'm|n(l,2) 

=  -s 4>n  sAtm]n 

(24) 

Cm|n(l.  3) 

=  c.(j>n  sA£mln 

(25) 

Cmi"(2. 1) 

=  s <pm  sA^m|n 

(26) 

Cm|n(  2,2) 

=  s<frm  s 4>n  cAf71'"  +  c 4>m  c <f>n 

(27) 

Cm  ln(2,3) 

=  -  s 4>m  c<f>n  cAC"|n  +  c(f>m  s<j>n 

(28) 

Cm|n(3,l) 

=  -c0n  sAC"1” 

(29) 

Cm|"(3, 2) 

=  -c<t>m  S(j)n  cAlm\n  +  s <j)m  c 4>n 

(30) 

Cm|n(3, 3) 

where  AHn  -  t1  -  im . 

=  c 4>m  c<t>n  cA£mln  +  s <T  s<pn 

(31) 
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Coordinate  Systems  | 

The  elements  of  Dm\n  are  given  by 

Dml”(l)  =  (Nn  +  nn)  c<f>n  sA^n  (32) 

Dm|n(2)  =  -(Nn  +  nn)  c <\>n  s <f>m  cAf”|n 

+  (Nn{l-e2)+nn)c<t>ms(j)n  +  Nme2C((>rns4>m  (33) 

Dm|n(3)  =  ( Nn  +  Hn)  c <t>n  c 4>m  cA£mln  -  (Nm  +  Um) 

+  {Nn(  1  -  e2)  +  Hn)  s 4>m  s<t>n  +  Nme2(s<f>m)2 


(34) 
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Coordinate  Systems  k 

The  deck  frame  is  colocated  with  the  ship’s  ENU  frame  with  the  x-axis  starboard,  the  y-axis 
forward,  and  the  z-axis  up.  This  sequence  of  plane  rotations  for  roll,  pitch,  and  heading  for 
platform  n  denoted  by  9™ .  6™ .  and  is  given  by 


c®n|n(i,  i)  =  c o%  +  senr  se;  s0?x 

(35) 

CVn |n(l,2)  =  -  C(9"  s0£  +  s(?"  S0;  c,0nh 

(36) 

CVn  l"(l,3)  =  s6?cff£ 

(37) 

cDn|n(2  i)  =  c0n  s0n 

(38) 

CVn\n(2. 2)  =  c6^c6l 

(39) 

CVnln(  2,3)  =  -s0" 

(40) 

CVn]n( 3, 1)  =  -  s 0"  c0l  +  s0”  s0£ 

(41) 

CVn\n{ 2,  3)  =  s0"  s(9,”  +  c0”  s6>”  c0" 

(42) 

C®n|n(3,3)  =  c0”  c0; 

(43) 
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Coordinate  Systems^ 


For  two  platforms  with  positions  V*  and  relative  to  Te  and  reference  coordinate  frames  Tn 
and  respectively,  the  position  of  target  T  relative  to  platform  n  is  denoted  by  position 
vector  1Zn  and  the  corresponding  position  vector  from  platform  m  is  given  in  deck  coordinates 
by 


where 


_  QDm\m'j^rn  _  £»X>m|n^n  jjT>m\n 


QDm  |n  _  £Vm\m£*m\n 

jjVm.\n  _  £tT>in\m  j-)m\n 


(44) 


(45) 

(46) 
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Coordinate  Systems~|| 


For  two  platforms  with  positions  V *  and  V ^  relative  to  Te  and  reference  coordinate  frames  Tn 
and  T™ ,  respectively,  the  position  of  target  T  relative  to  platform  n  is  denoted  by  position 
vector  7 Zn  and  the  corresponding  position  vector  from  platform  m  is  given  in  deck  coordinates 

by 


where 


(47) 

QVm\n  _ 

^jVm\m  Qm\n 

(48) 

jyDm\n 

^'Drn\m  jjrn\n 

(49) 
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The  Cartesian  coordinates  of  a  position  vector  in  Jrm  are  given  in  terms  of  the  spherical 
coordinates  by 

rm  cos  em  sin  bm 

r m  cos  em  cos  bm  (50) 

rm  sin  ern 

Because  the  deck  frame  is  colocated  with  the  local  geographic  frame,  the  sensors  may  measure 
the  range  rVn ,  bearing  bVn  and  elevation  eVn  in  deck  coordinates. 


Georgia  Tech  Research  Institute 


31 


Coordinate  Systems 


From  (50),  the  spherical  coordinates  are  given  in  terms  of  the  Cartesian  coordinates  by 
rm  =  \/  (xm)2  +  ( ym )2  +  ( zm )2 


■’(9-  * 


2  Vx 


'(I;).  I3-"1 1  >  \ym\ 


/(xm)2  +  (ym)2 1  >  | zm| 


V(^m)2  +  (ym)2/ 

;-ta„-(v^5+55!),  |,.|>|^v  +  (F)'ai 


The  radial  velocity  or  range  rate  is  given  by 


xmxm  +  ymym  +  zmzm 
y/(xm)2  +  {ym)2  +  (zm)2 


Target  Motion  Model  I 

The  dynamical  equation  that  is  commonly  used  to  represent  the  motion  of  the  target  relative 
to  the  platform  is  given  by 


Xfc+1  -  FkXk  +  Vk 


where 


Xk  —  target  state  vector  at  time  4 

Fk  =  linear  constraint  on  the  target  dynamics  at  time  tk 

Vk  —  white  noise  error  in  the  state  process  with  14  ~  Af( 0,  Qk) 


The  N(x.P)  denotes  the  Gaussian  distribution  with  mean  x  and  covariance  P. 


Target  Motion  Model 


The  state  vector  for  a  nearly  constant  velocity  motion  model  is  given  by 


Xk  = 


%k  Xk  yk  Uk  Zk  zk 


where  (xk-,yk,Zk)  represents  the  position  of  the  target  in  Cartesian  coordinates  at  time  tk  and 
(±k,Vk:  Zk)  represents  the  velocity  of  the  target.  The  dynamical  constraint  for  nearly  constant 
velocity  motion  is  given  by 


- 1 

?r 

02x2 

02x2 

where  Ak  — 

1 

6k 

02x2 

Ah 

02x2 

0 

1 

02x2 

02x2 

L 

5k  =  tk  —  tk- 1  denotes  the  sample  period,  tk  denotes  the  time  of  measurement  k,  and  02x2 
denotes  a  two-by-two  matrix  of  zeros. 


One  approach  for  filter  design  is  to  set  y/q5k  approximately  equal  to  the  maximum  acceleration 
of  the  target  under  track. 
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EKF  for  Radar  Measur^mentsJ 

The  measurement  equation 

is  given  by 

Z*  =  h,(Xk)  +  Wk 

(60) 

where 

Z*  = 

measurement  vector  at  time  tk 

= 

target  state  vector  at  time  tk 

wk  = 

white  measurement  error  with  Wk  ~  Af(0.Rk) 

(61) 
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where 


EKF  for Radar  Measurements^ 


= 


Tk 

bk 

ek 

Tk 


=  h(X  k)  +  Wh 


’  hr(Xk)  ' 

W rk 

hb(Xk) 

+ 

Wbk 

he(Xk) 

Wek 

_  hD(Xk )  _ 

wDk 

Tk 

bk 

ek 

Tk 

Wrk 

Wbk 

Wek 

WDk 


(62) 


measured  range  at  time  tk 

measured  bearing  at  time  tk 

measured  elevation  at  time  tk 

Doppler-derived  range  rate  at  time  tk 

error  in  range  measurement  with  wrk  ~  Af(0,  a^k) 

error  in  bearing  measurement  with  Wbk  ^  N(0.crlk) 

error  in  elevation  measurement  with  wek  ~ 

error  in  Doppler  range  rate  measurement  with  wok  ~ 
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EKF  for  Radar  Measurements 


hr(Xk) 

hb(Xk) 

he(Xk) 

hD(Xk) 


'xl  +vl  +  zl 


.-i (xk 


\yk' 

tan_1(  — ),  \xk\>\yk\ 

\xk J 


/xk  +  yk 


^k&k  d*  VkVk  “1“  Zk^k 

V4  +  yl T4 


xl+y\\>\zk\ 


1**1  >  \Vxt  +  vl\ 


EKF  for  Radar  Measurements  I 

The  EKF  for  the  target  state  estimation  with  nonlinear  measurements  is  given  by 

State  prediction: 

Xk\k-i  — 

Fk-iXk-i\k~i 

(67) 

Pk\k-i  = 

Fk-\Pk-i\k-iFk-i  +  Qk-i 

(68) 

State  update  with  the  measurement: 

Xk\ k 

=  Xk\k~i  +  KkZk 

(69) 

Pk\k 

=  [/  -  KkHk]Pk\k_x 

(70) 

Zk 

=  Zk-h(Xk  |jb_i) 

(71) 

Kk 

=  Pk^HtS,1 

(72) 

Sk 

=  HkPk\k-\Hk  +  Rk 

(73) 
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EKF  for  Radar  Measurements 


state  estimate  at  time  tk  given  measurements  through  tn 
filter  residuals  at  time  tk 

error  covariance  at  time  tk  given  measurements  through  t,j 

'  Hrk  ' 

[  dhr(Xk)  1 

dXk 

dhbjXk)  Hhk 

dXk 
dhejXk) 

dXk  Hek 

dhpjXk) 

dXk  J  xk=xk,„.. 


EKF  for  Radar  Measurements 


Xk\k-\  =  [  a 

the  rows  of  Hk  are  given  by 

Hrk  = 

rk\k- 

Vk\k 

Hbk  =  r/)2 

r,lk\, 

Hek  =  [ 

L  Tk\ 

HDk  =  [  HDk 


—  1  £fc|fc-l  Vk\k-\  Vk\k-\  %k\k-\  *k\k-\ 


1 

Vk\k_ 

—  o 

Zk\k 

—  0 

fk\k-i 

rk\k 

-1 

rk\k 

-l 

Vk\k-\ 

n  ^ 

0 

0  0 

rak\k~l 

r 

hk\k-l 

xk\k- 

i  Zk\k—i 

0  - 

Vk\k- 

-1  zk\k- 

r2 

f  k\k-i 

rhk\k-l 

r 2 

1  k\k- 

xrhk\k- 

HDk(  1) 

xk\k-l 

HDk(  3) 

Vk\k-\ 

rk\k-l 

HDk(5) 


Zfc|fc-i 

r’fcIA:-  1 


EKF  for  Radar  Measurements 


The  covariance  of  the  measurement  errors  in  the  platform  reference  frame  is  given  by 
tVrk  Hrk  ®rk  ^  ^ 

COV  Wbk  =  Hbk  MlQEi  o  a\k  0  QEsk)T Mk  [  HJk  H?  Hjk  )ss) 

.  Wek  J  [  Hek  °  0  a\k  _ 

where  a^k,  &bk  and  a\k  are  the  variances  of  the  measurements  in  the  sensor  reference  frame  and 

r  —  1  %k\k—\Zk\k—\  " 

-  Vk\k-l - ~T - 1 - 

rk\k-i  rnk\k-i 


Et  =  Ek 


+Lk 


with  J Ek  = 


Hk\k-\ 


yk\k-\zk\k-i 


%k\k~  1  L 

^fe|*:-l  r^k\k-l 


2fc|fc-l 

rk\k-l 


1  0  0  0  0  0 

QT  =  001000 

0  0  0  0  1  0 
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EKF  for  Radar  Measurements 


When  considering  the  measurements  in  the  u-plane  and  the  v-plane. 


Z*  = 


hr(Xk) 

hu(Xk) 

hv(Xk) 

hD(Xk) 


where 


uk  =  measured  sine  of  the  angle  of  the  target  in  the  ?t-plane  at  time  t.k 

vk  =  measured  sine  of  the  angle  of  the  target  in  the  v-plane  at  time  tk 

wuk  =  error  in  the  uk  measurement  with  wuk  ~  Af(0,  < 72uk) 

wvk  —  error  in  the  vk  measurement  with  wvk  ~  fif( 0,  cr2vk) 

hu(Xk)  =  .  ,  **., . == 


hv{Xk)  = 


'x 1  +  vl  +  4 


/xl  +  vl  +  zi 


EKF  for  Radar  Measurements 


Covariance  of  the  measurement  errors  in  the  platform  reference  frame  is  given  by 


Wrk  Hrk  &rk  ^  ^ 

COV  Wbk  =  Hbk  MlQUi  0  alk  0 
^ek  _  _  ffefc  _  ^  ^vk 

*  (QVi? M„[  H?k  H*  Hi 

where  a2,,  a2* .  and  rr2.  are  variances  of  the  measurements  in  the  antenna  aperture  and 


u£  =  uk 


uk  = 


I  Xk\k~  i  —  Xk\k-  \  4"  Lk 
xk\k-l 

—  -  rk  !*_! 

rk\k~\ 

Vk\k-\  ^ 

rk\k-i 

zk\k—l  xk\k-l 

—  - - 1 - rk\k 

zk\k-i 


Tk\k-\ 


rk\k-l 


Vk\k~i 

Zk\k-1 


Tk\k-1 
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^rGcessi^^ 

For  local  measurements  on  platform  m,  the  dynamical  motion  and  measurement  equations  are 
given  by 

x?  =  ^  +  vt-1  (94) 

Zk  =  hm(Xr  +  M$^Li^)  +  B?  +  WT  (95) 

where 

X ™  =  target  state  vector  in  ENU  on  platform  m  at  time  tk 

Z ™  —  measurement  vector  at  time  tk  on  platform  m 

F™  =  constraint  on  target  dynamics  at  time  tk  on  platform  m 

—  rotational  from  ENU  on  platform  m  at  tk- 1  to  ENU  on  platform  m  at  time  tk 

M-lkVrn  =  rotation  from  deck  coordinates  of  platform  m  at  tkto  ENU  on  platform  m  at  tk 
Lk^k  —  translation  from  ENU  on  platform  m  at  time  tk- i  to  ENU  on  platform  m  at  time  t 

LSk^Vrn  —  translation  in  deck  coordinates  from  ENU  on  platform  m 
at  tk  to  the  sensor  frame  on  platform  m  at  tk 
14™  1  =  error  in  the  state  process  for  platform  m  with  14™  1  rsj 
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Processing  Local  Measurements 


W™  =  measurement  errors  for  platform  m  with  W™  ~  A/r(0, R™) 
B™  =  residual  sensor  bias  for  platform  m  with  B ™  ^  A/*(0,  i?f  m) 


Processing  Local  Measurements 


'fcifc-1  U3x3 


where  denotes  the  rotation  of  (21)  from  ENU  on  platform  n  at  time  tj  to  ENU  on 
platform  m  at  time  tk-,  and 


1  0  0  0  0  0 
0  0  1  0  0  0 
0  0  0  0  1  0 
0  1  0  0  0  0 
0  0  0  1  0  0 
0  0  0  0  0  1 


<-i(3) 


where  Dk denotes  the  translation  of  (22)  from  ENU  on  platform  n  at  time  tj  to  ENU  on 
platform  m  at  time  tk . 


Processing  Local  Measurements! 

EKF  for  the  target  state  estimation  on  platform  m  with  nonlinear  local  measurements  is  given 

by 

State  prediction: 

\rm  _  A/f171  lm  Z?m  vm  I  rm!m  /no\ 

Ak\k-\  ~  Mk\k-l*k-lAk-l\k-l  +  Lk\k-1 

Pfik-i  =  +  Qk-x  (99) 

State  update  with  the  measurement: 


\rm  _  i  IS™  7™ 

Ak]k  —  Ak\k-l  +  ^k  Lk 

K  =  Zr~hm(X^k_1  +  M^mL^Vm) 

nrn  _  f  r  Tsm  rrmi  pm 

“k\k  —  U  —  Kk  ttk 

KJ?  =  Pfik-AWFlsr}-1 

^k  “  Mk  ^k\k-\\Hk  )  +  nk  +  rik 

Z ™  =  filter  residuals  on  platform  m  at  tk 

dhm(xk)  | 


(105) 


^Proxjessi^ 

For  local  measurements  in  deck  coordinates  on  platform  m,  the  measurement  equation  becomes 
zvm  _  hm(Mk™]mX?  +  LSk™]Vrn)  +  BJ?  +  W£l  (106) 

where 

—  rotational  from  ENU  on  platform  m  at  t k  to  deck  coordinates  on  platform  m  at  tk 
For  processing  the  measurements  in  deck  coordinates,  the  measurement  update  is  given  by 


y  m 

A k\k 

y  m  ,  Ts'Dm  /yT>m 

~  Afc|fc-1  +  AA:  Lk 

(107) 

rjrT>m 

ry'Drrt  t.  m  /  ■\.r'Dm\m.  y  rn  ■  j  Sm\T>m\ 

-  -  h  (Mk]k  +  Lklk  ) 

(108) 

-pm 

r k\k 

=  [/  -  K^H^M^]PSk-, 

(109) 

Kkm 

(110) 

iT>m\m 

k 

=  +  R? 

(111) 

where 


dh™{Xk) 

dXk 


X^M^X^+Ll^ 


Processing  Local  Measurements 


The  covariance  of  the  measurement  errors  in  ENU  for  platform  m  is  given  by 


0  0 


hs  Kr'Te^"""*  o  o 


0  0  si 


x  [  (HS)T  (flff)T  (H3)’ 


where  a^k.  v\k,  and  a\k  are  variances  of  the  measurements  in  the  sensor  reference  frame  and 


ntSmjm  _ 


\xkik_i=M*™'mXZlk_l+Lsk™'Vm 


with  Ek  defined  by  (86)  and  Q  defined  by  (87). 


Processing  Local  Measurements! 

For  angular  measurements  in  terms  of  u  and  v  as  in  (88).  the  measurement  covariance  matrix 
for  range,  bearing,  and  elevation  in  ENU  is  given  by 

’  w™k  1  H?k  1  d2rk  0  o' 

cov  „s  =  HS,  nCl”)TWfH”  0  clk  0 

.  •»;;  II//,';  loo  _ 


x  [  (ffs)1  (hs)t  (ff;;)T 


where  and  o\k  are  variances  of  the  measurements  in  the  sensor  reference  frame  and 


j^jStTi\tTL  _ 

k  k  v  _U5mlm  vm  ,  r-Smj'Dm 


with  Uk  defined  by  (93) 


Processing  Remote  Measurements! 

For  processing  remote  measurements  from  ENU  on  platform  n  in  ENU  on  platform  m,  the 
dynamical  motion  and  measurement  equations  are  given  by 


K  =  ''C”  f"-  ,  -Y” .  +  .  +  v~- 1 

Z?  =  ft”[Mr,rxr+z;|r+M^"ifjn+G;l'n  +  c;  +  H'f 


where 


Qn\m  _  resj(juai  gridlock  bias  for  ENU  on  platform  n  relative 
to  ENU  on  platform  m  with  ~  J\f( 0,  Rk  n^m) 

The  exact  form  of  and  R^nlrn  depends  on  the  method  of  gridlock  (e.g.,  relative  or 
absolute). 
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Processing 

EKF  for  target  state  estimation  in  ENU  on  platform  m  with  with  nonlinear  remote 
measurements  in  ENU  on  platform  n  is  given  by 

State  prediction: 

y  m  _  Z?m  Ym  I  T  mIm 

Ak\k~X  -  Mk\k-lPk~l^k~l\k~l  + 

nm  _  7i /fm lm  T>rn  (  A/fmlm  Z?m  ,  yem 

State  update  with  the  measurement: 


ym 

Afc|fc 

ym  i  tv'TZ 1 7Ti  yn|m 

~  Afc|fc-1  + -^Jfe 

^njm 

=  -  hn  (Mkj™  Xk\k_1 

pro 

11 

1 

j?3 

*•  s 

£ 

^n|m 

r<n]m 

Trnjm 

dhn(Xk ) 
dXk 

* n\T>n  rSnlDn^ 
'Lk\k  ^k\k  ) 


,+^!:+<nC" 
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Processin^Remot^MeasurBmervt^r^ecl^i^m^| 

For  remote  measurements  in  deck  coordinates  on  platform  n,  the  measurement  equation 
becomes 

zf"  -  /."(Mj”i”xr+Mjr'Xi'r+c>'I’")+s2+M»l"G?m+w'‘  (i2?) 

<i28> 

For  processing  measurements  in  deck  coordinates  on  platform  n,  the  measurement  update  is 
given  by 

X%k  =  X%k_r  +  K”nlmZkn]m  (129) 

Z?n\m  =  Zfn  _  jqn  +  +  Lf”115”)  (130) 


z?"  -  k'w^xfr., + + ifs1®") 


=  [7  —  Kk  Hk 


Vn\m  TjVn\m  ,  fX>n|mn  pm 


'jUfc  fcj/c  — 1 


Kvknxm  =  pfcTfc_1(^n|m<r|m)T[^n|m]-1  (i32) 

5®n,m  =  pr|mM?r|m^-i(^r|mMS|m)T + Pfcn^r|nj?fn|m(MS|n)r + ^33) 

jj-Pnlm  _  9h  {Xk)  |  ('1^4^ 


+  M®"l”L"!r  +  Lfl"10" 


fc|/c  *1*  ^  k\k 
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Processing  Time-Delayed,  Out-of-Sequence  Measurements 


Processing  the  time-delayed,  out-of-sequence  measurement  Z^-i  is  given  by 
State  update  with  the  measurement: 

Xk\k  =  Xk]k\k-\  +  Kk-\\Zk-i  -  h{F^(Xk\k\k_i  -  QkHl S^'Zk)))  (135) 

Pk\k  =  [I-Kk^Hk-XF^}Pkm_,  +  Kk^Hk„xF^Qk{I-KkHk)T  (136) 

Kk- 1  =  [Pk\k\k-i-{I-KkHk)Qk]{Hk-XF^)TS^\  (137) 

Sk-i  =  Hk-1F£1[Pmk_l-(I-KkHk)Qk-Qk(I-KkHk)T  +  Qk 

-  QkHlS^HkQk){Hk^F^)T  +  /?*_!  (138) 


where 


Kk  = 

Kalman  gain  at  tk 

zk  = 

filter  residuals  at  tk 

5,  = 

covariance  of  the  residuals  at  tk 

Hk  = 

linearized  output  matrix  at  tk 

Hk- 1  = 

linearized  output  matrix  at  tk- 1 

Note  that  <5*.  =  tk  ~  is  positive  in  the  definition  of  Fk  and  Qk  as  specified  by  (58)  and  (59). 
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Interacting  Multiple  Model  (IMM)  Estimator 


The  dynamical  motion  and  measurement  equations  for  system  with  Markovian  switching 
coefficients  and  nonlinear  measurements  arc  given  by 

Xk+\  =  Fk(Qk+\ )Xk  +  Vk(0k+\)  (139) 

Z*  =  h(0k,Xk)  +  Wk  (140) 

where  9k  is  a  finite-state  Markov  chain  taking  values  in  {1, . . . ,  N]  according  to  the  transition 
probabilities  Pij  of  switching  from  mode  i  to  mode  j,  and 

E[Vk(8k+i  -  i)V?(8k+ 1  =  i)}  =  Qk(6k+ 1  =  0- 
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Concluding  Remarks^ 


Multiplatform-Multisensor  Data  Association 
Track  Management 

-  Track  Number  Management 

-  Split  or  Diveraging  Tracks 

-  Track  Swaps 

-  Track  Deletion 

-  Deleting  Redundant  Tracks 
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